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ABSTIV^CT 


The  Walsh  Transform;  is  invsstigatsd  for  its  usefulness  in 


deinterleaving  the 


inte^l 

ESMjr« 


eaved  pulse  stream  presented  to  a 


preprocessor  by  an  E^i^'receiver .  After  background  chapters  on 
a  typical  ESM  system  and  the  theory  and  characteristics  of  the 
Walsh  Functionsi  a  number  representation  of  the  pulse  stream 
is  described.  Fast  Walsh  Transforms  and  Power  Spectral 
Densities  of  the  pulse  representations  are  computed  and 
analyzed  for  features  that  could  be  used  to  recognize 
individual  pulse  trains  in  the  interleaved  representation. 
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IhlTRODUCTION 


On  a  June  night  in  1940,  a  lone  ANSON  aircraft  -flew  a 
solitary  path  through  the  sKies  of  East  Anglia,  England.  The 
wireless  operator  was  listening  for  traces  of  signals 
propagaf'ng  through  the  dark  space  around  the  aircraft.  Near 
an  important  military  target,  a  signal  came  through  loud  and 
clear  in  his  earphones.  First,  a  series  of  dots  at  the  rate 
of  sixty  per  minute.  Then,  as  the  ANSON  droned  onward,  the 
dots  disappeared  into  a  steady  tone.  A  little  later,  the 
steady  tone  broke  up  into  a  series  of  dashes,  at  the  rate  of 
sixty  per  minute. 

This  signal,  known  as  a  Lorentz  beam,  was  transmitted  from 
Germany  and  it  posed  a  considerable  threat  to  the  high  value 
military  target.  The  signal  was  really  composed  of  two  beams, 
one  beam  of  dots,  and  one  beam  of  dashes,  shown  in  Figure  1.1. 
By  using  the  reception  of  the  dots  and  dashes  as  an  indication 
of  whether  they  were  on  course,  German  bombers  could  navigate 
to  their  targets.  The  phasing  of  the  dots  and  dashes  was  such 
that  they  combined  into  one  steady  tone  in  the  overlap. 

Receipt  of  the  steady  tone  meant  the  aircraft  was  on  the 
course  laid  down  by  the  transmitting  site  and  receipt  of  dots 
or  dashes  meant  the  bomber  was  right  or  left  of  the  beam 
overlap.  With  that  information,  the  bomber  crews  could  easily 
get  back  on  course  and  perform  their  mission.  Cl] 


The  beeine  were  ultimately  countered  by  the  British  using  a 
variety  of  methods,  but  the  invaluable  intelligence  gained  by 
the  bomber  over  East  Anglia  vuas  only  the  beginning  of  a 

wartime  effort  that  was  to  save  thousands  of  lives,  both 
civilian  and  military.  That  brave  and  skillful  act  has  been 
repeated  hundreds  of  times  since,  in  a  'police  action”  and  an 
undeclared  war,  but  the  priceless  results  have  been  the  same, 
that  of  saving  lives. 

Key-  Stei^j  note  zone 

Trjnsmitter 


Dots  tfishd  Interlock 


so  that  Mrhcre  beams  overlap 
a  steady  note  is  hear^ 


THE  LOI^ENZ  BEAM 


Figure  1.1.  The  Lorentz  Beam 

Since  the  arrival  of  the  first  electronic  transmitter  and 
receiver,  man  has  been  listening  to  his  own  radiations  for 
coiranuni  cat  ions,  pleasure,  and  information.  Norld  Mar  II, 
however,  brought  a  new  maturity  to  listening  to  his  own 
radiations.  Now  he  was  listening  to  an  opponent's 
transmissions  to  derive  and  use  information  about  them. 
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About  1969 f  the  term  Electronic  Support  Measures  <£SM)  was 
coined  C2:  p.S71  and  now  reflects  the  important  business  cv 
obtaining  and  using  information  gained  on  the  electromagnetic 
radiations  of  another.  The  alert  operator  aboard  the  lone 
British  aircraft  used  his  early  and  undefined  ESM  receiver  to 
obtain  Knowledge  that  was  used  against  the  German 
transmi tters. 

A.  THE  ELECTRONIC  l^RFARE  <EW)  ENVIROMIENT 

Of  all  the  sensors  that  an  enemy  might  use,  the  most 
important,  at  least  as  an  electronic  threat,  is  radar.  This 
sensor  is  almost  always  the  long  and  short  range  eyes  and  ears 
of  the  opposing  force,  and  denial  of  its  use  gives  a  strategic 
and  tactical  advantage  to  the  person  capable  of  doing  so.  I2t 
p.71 

Consider  the  case  of  a  penetrating  bomber.  In  order  to 
Jam,  deceive,  or  otherwise  negate  the  effectiveness  of  the 
radars  of  the  defense  system,  emissions  must  be  detected, 
analyzed,  and  classified  as  to  their  threat.  In  this 
endeavor,  the  ESM  receivers  being  used  will  not  lacK  for 
electromagnetic  sources  to  perform  their  analyses.  Early 
warning,  target  detection,  acquisition,  and  tracking,  or  even 
missile  and  interceptor  guidance  and  control  radars  will  paint 
the  aircraft  with  an  invisible  illumination. 

These  electromagnetic  energies  will  be  of  many 
frequencies,  antenna  scans,  pulse  shapes,  and  modulations,  and 
even  though  radar  systems  are  well  understood  and  sometimes 
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easily  countered,  the  sheer  numbers  and  types  o-f  received 
signals  mill  tax  the  capabilities  o-f  the  £SM  system. 

As  the  bomber  nears  the  target,  the  signal  densities  Mill 
soar  because  the  de-fense  is  toughest  close  to  the  target.  For 
instance,  it  has  been  reported  that  over  168  SAM  sites  were 
once  located  on  the  Egyptian  side  o-f  the  Suez  Canal.  As  each 
battery  transmits  several  signals,  easily  308-508  signals 
could  be  -found  -from  this  one  class  o-f  weapon.  Mhen 
communication  and  radio  location  signals  are  added,  the  signal 
density  problem  is  even  -further  complicated.  12:  pp. 27-331 

Quite  easily,  the  combat  environment  that  the  bomber  is 
in  will  reach  566,980  pulses  per  second.  C3:  p.S41 

Figure  1.2  141  gives  an  indication  o-f  the  pulse  densities 
the  bomber  might  encounter.  The  pulse  density  is  dependent 
upon  both  system  sensitivity  and  platform  altitude  but  is  not 
uni-formly  dependent  on  these  parameters.  At  low  sensitivities 
<-20  to  -70  dBm)  pulse  density  varies  mostly  with  sensitivity 
and  not  with  altitude.  At  high  sensitivity  0-100  dBm)  pulse 
density  is  almost  totally  altitude  dependent. 

The  identi-f ication  o-f  various  emitters  and  the  evaluation 
o-f  the  possible  threats  they  pose  to  the  aircraft  requires  the 
sorting  of  these  pulses.  This  separating  or  deinterleaving 
of  signals  is  a  difficult  problem  that  is  discussed  in  Chapter 
II  . 
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B.  ELECTRONIC  WARFARE  SUPPORT  MEASURES  (ESM> 


Th«  broad  term  Eloctronic  War-fare  (EW>  ,  or  sometimes 
referred  to  as  Electromagnetic  Warfare  (EMW> ,  is  used  to 
denote  military  actions  that  involve  the  friendly  exploitation 
of  the  electromagnetic  spectrum  and  the  prevention  of  hostile 
use  of  the  el ectromagetic  spectrum.  It  should  be  realized 
that  the  electromagnetic  spectrum  covers  from  zero  frequency 
to  the  infinite  frequency,  which  includes  the  optical, 
infrared,  and  laser  frequencies.  ESI  Coverage  of  these 
additions  to  the  RF  spectrum  adds  considerably  to  the  problem. 

£W  may  be  divided  into  three  major  divisions: 

1.  Electronic  Warfare  Support  Measures  <ESM) 

2.  Electronic  Counter  Measures  (ECM> 

3.  Electronic  Counter  Counter  Measures  <ECCM> 

Elaccronic  Support  Measures,  the  usual  shortened  form  of 
Electronic  Warfare  Support  Measures,  will  be  defined  shortly. 
ECM,  quite  briefly,  is  the  generation  of  intentional 
electronic  interference  between  electronic  systems  for  a 
useful  purpose.  ECCM  ,  also  quite  briefly,  is  action  taken  to 
counter  the  detrimental  effects  of  the  opponent's  ECM. 
p.  13 

This  section  will  deal  only  with  the  ESM  arena. 

1.  A  Definition  of  ESM 

Consider  a  specific  definition  of  ESM.  ESM  is: 

That  division  of  Electronics  Warfare  involving  actions  to 
search  for,  intercept,  locate,  record,  and  analyze  radiated 
electromagnetic  energy  for  the  purpose  of  exploiting  such 
radiations  in  support  of  military  operations.  C7:  p.6] 


This  ci»+inition  carries  an  underlying  statement:  ESM 
is  the  source  o-f  in-formation  required  to  carry  out  the  other 
Electronic  Mar-fare  divisions  as  well  as  threat  detection, 
warning,  avoidance,  target  acquisition,  and  homing.  C4: 
pp.3~d]  In  addition,  by  examining  the  individual  actions 
speci-fied  in  the  de-finition  we  can  gain  -further  insight  into 
the  ESM  arena. 

The  search  -for  and  interception  of  electromagnetic 
radiations  is  simply  a  basic  part  of  the  ESM  mission.  The 
intercepted  signal  must  be  analyzed  and  a  determination  made 
of  its  character,  parameters,  and  location.  Locating  the 
emitter  is  usually  a  function  for  direction  finding  antennas. 
Recording  of  the  signal,  certainly  the  parameters,  is  usually 
done  simultaneously  with  the  analysis.  On  some  missions  a 
recording  is  made  of  the  signal  for  postflight  analysis  with  a 
mimimum  of  information  being  displayed  to  the  operator  in  real 
time.  Generally  it  is  desired  to  produce  this  information 
from  an  ESM  systim  with  an  inherent  automatic  processing 
capability  that  detects,  classifies,  and  flashes  a  warning  if 
need  be.  12 t  p.36] 

The  information  that  the  ESM  mission  provides  ranges 
from  that  collected  with  regularly  scheduled  and  dedicated 
collection  excursions  to  debriefs  of  crews  after  completion  of 
a  routine  patrol  mission.  To  a  tactical  commander,  this 
information  is  more  extensive  and  useful  than  plain 
intelligence.  It  provides  the  precise  state  of  the  electronic 


including  technical  characteristics  and  emitter 
location,  and  becomes  part  o-f  mission  pre-flight  brie-fingp. 

In  obvious  -fact,  radars  not  transmitting  are  radars 
that  can't  be  intercepted  and  analyzed.  Thus,  i-f  there  are  no 
transmissions  from  specific  desired  radars,  or  any  at  all, 
tactics  that  might  provoke  them  into  radiating  are  viable  in 
the  ESM  scenario.  Indeed,  this  action  is  included  in  the 
basic  rules  for  Electronic  Reconnaissance.  t2:  p.65I 
2.  The  Role  of  ESM  in  Electronic  Warfare 

In  general,  ESM  is  passive  electronic  warfare,  ECM  is 
Active  electronic  warfare,  and  ECCM  may  be  either.  These 
distinctions  may  not  be  clear  cut,  however.  ESM  may  involve 
the  radiation  of  a  signal  to  determine  the  characteristics  or 
electronic  reaction  of  enemy  equipment.  ECM  may  involve  the 
passive  reception  of  enemy  signals  in  order  to  decide  which 
one  to  counter  C2:  p.3].  The  latter  could  be  the  use  of  a 
tactical  ESM  system  integrated  with  an  ECM  system  for  platform 
protection.  The  relationships  between  ESM,  ECM,  and  ECCM  are 
shown  in  Figure  1.3.  C2:  p.2] 
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Figure  1.3.  The  Interactions  of  EW 


The  previous  consideraition  o-f  the  penetrating  bomber 
can  be  continued  in  the  context  o-f  ESM's  role  in  electronic 
warfare. 

The  bomber  will  almost  surely  face  an  integrated 
defense  network.  It  is  absolutely  essential  to  Know 
beforehand  the  technical  characteristics  and  location  of  the 
electronic  defense  that  will  be  used  against  the  bomber.  Some 
Kind  of  electronic  reconnaissance  information  is  needed  in 
order  to  maximize  the  effectiveness  of  ECM  and  ECCM  equipment. 
Electronic  intelligence  <ELINT)  and  ESM  are  our  sources  of 
information . 

The  value  of  the  knowledge  gained  from  electronic 
reconnaissance  of  the  defense  system  is  proportional  to  its 
currency.  The  use  of  information  from  electronic  intelligence 
which  was  gained  from  many  sources  and  ov^er  a  longer  period  of 
time  is  important,  but  the  immediate  needs  of  the  tactical 
commander  are  met  using  ESM.  This  information  is  collected  in 
the  hours  before  takeoff,  undergoes  minimal  analysis,  and  is 
used  for  mission  planning.  C2:  p.373  Thus,  ESM  is  electronic 
reconnaissance  to  determine  the  present  state  of  the  defense 
system. 

ELINT,  which  has  been  subjected  to  extensive  analysis, 
has  been  used  to  design  and  develop  the  ECM  suite  of  the 
bomber.  But  the  tactical  commander  still  uses  ESM  to  plan  and 
accomplish  the  mission.  Though  they  may  share  a  common 
collection  platform,  the  difference  between  them  can  be  seen 
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in  thoir  op«raitiona1  control  ,  *0(1  use,  and  time  or  duration  of 
collection.  ESM  is  used  by  the  tactical  commander  today,  and 
ELINT  goes  to  the  rear  echelon  for  further  study.  C2:  p.743 
Figure  1.4  C83  shows  the  operational  role  of  ESti  in 
electronic  warfare.  Note  the  central  position  of  ESM  as  i t 
"feeds"  both  ECM  and  ECCM  and  operational  readiness. 


Figure  1.4.  Operational  Relationships  in  EN 


C .  SUMMARY 

The  previous  sections  only  touch  upon  the  vast  nature?  of 
EW  and  ESM.  They  provide  a  general  groundwork  for 
understanding  the  purpose  of  using  ESM  and  its  role  in  the 
electromagnetic  conflict. 

This  conflict  is  extremely  important.  New  and  modern 
weapon  systems  depe^nd  upon  victory  in  the  electromagnetic 
conflict  as  a  prerequisite  for  victory  in  battle.  [2:  p. 1] 

Underscoring  this  premise  are  new  policies  and  procedures 
by  the  Defense  Department  which  will  accelerate  development 
and  introduction  of  electronic  warfare  systems.  The  rapid 
escalation  and  introduction  of  new  Soviet  radars  places 
current  EW  sytems  at  a  disadvantage,  especially  in  countering 

i 

fire  control  and  missile  guidance  radars.  C93 

The  future  probably  holds  a  greater  integration  of  ESM 
<certain1y  tactical  ESM)  and  ECM  systems,  with  design  for 
military  aircraft  leading  the  way.  Already  systems  that 
combine  ECM  capability  with  passive  warning  equipment  for 
aircraft  are  being  sought  by  many  NATO  nations,  such  as  the 
British  Royal  Air  Force  C101,  the  Belgian  Air  Force  CHI,  as 
well  as  the  United  States  Air  Force  C 123 .  In  addition,  and 
with  considerable  impact  on  system  design,  a  trend  toward 
carrying  these  systems  internally  rather  than  in  pods  hung 
underneath  the  aircraft  seems  to  be  evident  C 10  and  123. 

In  the  hardware  area,  the  use  of  embedded  computers  and 
microprocessors  with  high  use  of  MHSIC  technology  will  become 
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prevalent,  along  with  all  o-f  the  ramifications  of  computer 
use.  The  sophisticated  signal  processing  needs  of  EH  will  have 
to  be  met,  with  verifiable,  reliable,  maintainable  software, 
and  the  development  of  high  capacity,  non-volatile  data 
storage.  C 13] 

Upcomin  and  foreseeable  developments  are  EH  systems  for 
space  applications  C2:  pp. 171-178].  All  of  the  unique 
environmental  conditions  of  space  will  generate  new  and 
different  challenges  for  ECM  and  ECCM.  These  must  be  met  in 
order  for  manned  space  stations  and  satellites  to  survive. 

D.  SCOPE  OF  THESIS 

( 

After  these  introdMctory  sections  Chapter  II  examines  some 
of  the  system  design  requirements  jf  an  ESM  system.  ~he 
various  components  are  looked  at  in  some  detail.  Signal 
processing  is  addressed,  beginning  with  basic  parameter 
measurement,  and  ending  with  a  discussion  of  our  basic 
problem,  that  of  deinterleaving  received  pulse  trains. 

The  attack  on  the  problem  begins  in  Chapter  III,  with  a 
discussion  of  the  theory  of  Walsh  and  Rademacher  Functions. 
Application  of  this  theory  is  carried  out  in  Chapter  IV,  where 
the  Walsh  transform  is  used  on  the  deinterleaving  problem. 
Finally,  conclusions  are  developed,  and  listings  of  the 
programt^  used  are  given  in  order  that  follow  on  work  can  be 
done . 

The  attack  on  the  problem  using  the  Walsh  Transform  will 
begin  with  basic  examples  and  proceed  with  an  application  of 


the  transform  directly  to  an  interleaved  pulse  train.  The 
goal  is  to  determine  whether  the  Walsh  Transform  outputs  a 
unique  feature  that  can  be  used  in  a  deinterleaving  algorithm. 

The  radar  pulses  can  be  represented  by  time  of  arrival 
<TOA)  tags,  which  are  output  by  the  ESM  receiver  to  the 
preprocessor.  A  TOA  tag  can  be  thought  of  as  merely  a  "marK" 
in  time  that  indicates  the  receipt  of  a  pulse.  The  stream  of 
incoming  pulses  forms  a  string  of  TOA  TAG'S.  The  Transform  is 
applied  to  the  TOA  tag  string. 

The  Walsh  Transform  is  easily  applied  to  binary  or  two 
level  functions,-  and  some  promise  is  held  for  its  use  in  a 
deinterleaving  algorithm  if  a  PRI/PRF  recognition  feature 
exists  in  the  transform. 

With  the  speed  and  ease  of  computation  of  the  Walsh 
Transform,  a  processor  based  on  these  functions  should  be  fast 
and  simple.  Consider  the  programs  used  to  compute  the  Walsh 
transforms  in  this  thesis.  They  were  originally  written  in 
FORTRAN,  and  were  easily  converted  to  BASIC  for  use  in  this 
thesis. 

This  simple,  easy  to  use  language  normally  brings  a 
response  that  would  question  the  use  of  such  an  "elementary* 
language.  That  response  is  misinformed. 

BASIC  is  easy  to  learn,  and  easy  to  use.  It  performs  well 
when  compared  to  other  structured  languages,  and  has  excellent 
mathematical  and  scientific  processing  capability.  The  use  of 
the  language  certainly  has  not  hindered  the  growth  of  the 
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computing  power  Available  in  a  microcomputer.  The  author^s 
personal  computer,  the  IBM  Personal  Computer,  with  6^K  o+  RAM, 
has  the  same  amount  o-f  memory  that  a  -frontline  U.S.  Navy  P3-C 
Orion  ASM  aircra-ft  uses  to  process  seven  tactical  work 
stations,  plus  handling  all  the  navigation,  displays,  and 
other  systems. 

Throughout  this  thesis,  a  microcomputer  was  used  to 
compute  all  the  results  given,  and  only  the  lack  o-f  a  thesis 
quality  graphing  device  prevented  the  author  -from  doing  all 
work  on  his  own  personal  computer.  Don't  underestimate  a 
microcomputer's  capabilities.  The  day  is  coming  when 
main-frame-like  computing  power  will  be  available  in  each 
o-f -f ice  as  close  as  the  desk  or  table. 


n.  ELECTRONIC  SUPPORT  MEASURES  SYSTEMS 


Th«  ESM  system  used  on  the  modern  bet  tie-field  must  be  able 
to  give  rapid  and  accurate  assessments  o-f  the  complex  RF 
environment.  Detection  and  recognition  of  hostile  radar 
transmissions  must  taKe  only  a  few  seconds  because  to  do  less 
invit'^s  destruction  or  the  loss  of  the  ability  to  deal  a 
decisive  blow  to  the  enemy.  C3:  p.54] 

The  design  of  ESM  systems,  in  view  of  their  importance, 
would  seem  to  follow  a  rational  and  orderly  path  from  concept 
to  development  to  operational  capability. 

This  usually  isn-'t  the  case,  unfortunately,  as  the 
reactive  nature  of  EM  will  impede  an  orderly  flow  or 
development.  Mi  thou t  perfect  knowledge  of  an  enemy's  radar 
and  EM  capabilities,  surprises  occur  which  tend  to  drive  the 
normal  development  into  a  "reaction*  driven  process.  The 
normal  process  also  tends  to  be  delayed  by  a  desire  to  wait 
until  a  precise  definition  of  the  threat  is  in  hand.  By  the 
time  it  is,  a  reaction  process  is  needed  to  counter  the 
threat.  C2t  p.93] 

As  referred  to  in  the  Introduction,  page  25,  this  design 
process  is  undergoing  changes  to  improve  the  response  and 
lessen  the  reaction  design.  Changes  in  the  design  process 
will  not  change  the  design  requirements. 


The  ESM  system  is  dlssignsd  to  intercept  many  different 
el ectrcxnsgnetic  signals.  Ideally,  the  system  should  be  able 
to  accomplish  as  quickly  as  possible  the  following  tasks: 

1.  Intercept  a  transmitted  signal  at  any  frequency 

2.  Determine  the  type  of  modulation  in  the  signal 

3.  Identify  the  usable  intelligence  carried  in  the  signal 

4.  Measure  the  direction  of  arrival  of  the  waveform  for 
calculations  that  locate  the  transmitter. 

9.  Process  and  record  the  signal  characteristics 
■  6.  Display  to  the  operator  and  transfer  to  the  computer  for 
decision  making 

In  short,  the  system  must  gather,  process,  and  display  all 
signals  of  interest.  C2:  p.l91] 

Meeting  these  requirements  requires  a  lot  of  attention  and 
hardware.  No  single  tuner  can  cover  the  area  of  interest 
(usually  900  MHz  to  30  GHz),  so  many  tuners,  antennas,  »nd 
other  system  units  are  required. 

Mith  this  equipment,  an  CSM  system  approaches  the  problem 
with  these  basic  intents: 

1.  Intercept  the  signal.  Conversion  to  a  usable  form  must 
now  take  place  immediately.  An  alarm  should  be  given  to 
indicate  a  signal  has  been  intercepted.  Some  analysis 
or  sorting  must  be  done  immediately  to  determine  if  the 
signal  is  of  high  threat,  since  signal  frequency  and 
modulation  usually  correlate  with  degree  of  threat.  12: 
p. 

2.  Analyze  the  signal  for  further  information.  Once  it's 
decided  there  isn^t  an  immediate  threat,  other 
parameters  of  the  signal  can  be  determined  for  recording 
and  later  decision  making. 

Each  of  these  approaches  has  more  or  less  importance 
depending  upon  the  platforms.  An  airborne  platform  has  less 
space  to  accommodate  equipment  for  reception  and  analysis,  but 
is  usually  exposed  to  more  immediate  threats.  Therefore,  an 


airborn*  plAt-form  dcvot««  mor«  »tt«ntion  to  intorceptlon  and 
iinn«diate  analysis  than  ths  ground  plat-form,  Mhich  probably 
is  mors  concsrnsd  with  a  full  analysis  and  rscording. 
Hopsfully,  this  sntirs  thssis  has  ths  flavor  of  applying  to  an 
airborns  platform. 

An  sxamination  of  a  typical  airborns  ESM  systsm  Mill  bsgin 
in  ths  nsxt  ssction. 

A.  COMPONENT  SYSTEMS 

Figurs  2.1  C4]  shows  a  typical  ESM  systsm.  Of  particular 
intsrsst  and  discussion  in  this  chaptsr  will  bs  ths  rscsiving 
systsm,  ths  prsprocsssor ,  and  ths  procsssor.  Ths  Data  Filss 
<Activs  Emittsr  and  Emittsr  Paramstsrs) ,  ars  ussd  by  ths 
systsm  as  rscords  for  currsnt  signals  bsing  rscsivsd  <Activs 
Emittsr)  and  comparison  of  paramstsrs  from  procssssd  signals 
with  prsviously  sntsrsd  <a  priori)  paramstsrs  of  signals  of 
intsrsst  or  sxpsctsd  signals  <Emittsr  Parsmstsrs) .  Elsctronic 
Countsr*-f1sasurss  <EGM)  handoff  and  squipmsnt  is  not  rsally  a 
part  of  ths  ESM  systsm,  but  is  includsd  as  a  function  that 
could  bs  avail abis  for  manual/automatic  jamming  or  dscsption 
of  high  thrsat  signals. 

1 ■  Ths  Rscsiving  Systsm 

tnssparabis  in  thsir  functions,  ths  antsnnas  and  ths 
rscsivsr  form  ths  intsrcsption  and  dstsr.tion  systsm  that  is  at 
ths  forsfront  of  ths  sntirs  procsss.  This  ssction  will  covsr 
a  fsw  of  ths  important  aspscts  of  ths  antsnnas  and  rscsivsr 
systems. 
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Figure  2.1.  An  ESM  System 


RECORDERS 


a.  Anti»nnA« 


Th«  «n«rgy  r«c«iv*d  by  «ny  platform  is  a  very 
smell  'frection  o-f  thet  transmitted.  Consider  Teble  2.1  C14t 
p.9'-33].  Notice  the  numbers  in  the  received  power  column  ere 
megnitudes  below  one  milliwett  in  el most  ell  instances. 


Teble  2.1. 

Received  Power 

Exempl es. 
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30 

9 
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18 

1 

38 

9 
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18 

298  K 

38 

18 

-24 

288 

-98 

46 

188  K 

33 

10 
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High  sensitivity  in  the  system  is  needed  to  use 
this  power.  Good  enter.ne  design  is  the  sterting  point  for 
reception  of  the  signels  thet  ere  desired  to  be  enelyzed*  end 
good  design  is  very  dependent  upon  the  epplicetion. 

The  besic  equetion  thet  is  fundementel  to  the 
emount  of  power  eveileble  for  use  ist 

Pp  -  PiG^X*/<<4lO  *R*>  <2.l> 

where  the  veriebles  end  units  ere  es  indiceted  belowt 


Pf.  «  r*c«iv»d  signal  powsr  in  watts 
•  transmittsd  powsr  in  watts 
6^  •  s-ffsctivs  transmitting  antsnna  gain 
Op  «  s-f-fsctivs  rscsiving  antsnna  gain 
X  >■  wavslsngth  <sains  units  as  R) 

R  ssparation  bstwssn  transmittsr  and  rscsivsr 

In  this  ags  o-f  *i-f  you  can  sss  it,  you  can  hit  it* 
accuracy  with  guidsd  missilss,  it  is  not  rscommsndsd  that 
shortsp  rangss  bs  ussd  to  improvs  ths  amount  of  powsr  bsing 
rscsivsd.  Thus,  ths  only  othsr  avail abis  factor  undsr  ths 
dssignsr^s  control  <bssidss  rscsivsr  dssign)  is  ths  rscsiving 
antsnna  gain. 

Our  antsnna  gain  is  affsctsd  by  ssvsral  factors, 

G  ■  4ltAP/X2  <2-2> 

whsrs  A  and  P  ars  dsfinsd 

A  *■  capturs  arsa 
P  *  antsnna  sfficisncy 

Ths  capture  arsa  is  ths  physical  sizs  of  ths 
antsnna.  The  sfficisncy  in  an  inhsrsnt  propsrty  and  a 
function  of  ths  typs  of  antsnna  <usua11y  bstwssn  8.9  and 
8.6)  . 

Ths  bsamwidth  of  "hs  pattsrn  that  has  ths 
dssignsT's  intsrsst  also  has  an  sffsct  on  ths  antsnna  gain 
C 193 .  It  is  an  invsrss  rslationship  t i4i  p, 29-31  that 
dspsnds  upon  ths  currsnt  distribution  across  ths  apsrturs. 

A  typical  rsf lector  antsnna,  for  instance,  has 
0  %  28,88e/<4B6B>  <2. 3) 

whsrs  80  and  #b  are  ths  half-power  bsamwidths  in  degress 
measured  in  ths  principal  planes  of  the  pattern. 


Th*  G*in  G  is  ths  powsr  gsin,  snd  should  be  used 
in  rsdsr  equations  because  it  includes  losses  introduced  by 
the  entenns  <Eq.  2.2).  Besring  resolution  and  coverage 
functions  for  direction  finding  would  lead  to  a  consideration 
of  the  directive  gain.  The  directive  gain  is  more  descriptive 
of  the  antenna  pattern. 

A  useful  expression  for  directional  gain  is 
Gq  •  4</(esSB>  (2.4) 

where  and  are  the  half -power  beamwiths  in  radians. 

C  13] 

Gain  must  be  considered  in  terms  of  signal 
reception.  Two  of  the  system‘'s  tasks  are  to  receive  the 
signals  which  may  come  from  any  direction  (omni 
considerations)  and  locate  the  emitter  (DF  considerations). 

The  two  tasks  are  not  necessarily  complementary,  so  a 
trade-off  in  system  performance  can  be  accepted  or  2  antennas 
can  be  used,  one  for  each  mission. 

For  an  airborne  platform,  weight,  space,  and 
structural  requirements  may  limit  the  size  of  an  ESM  antenna, 
but  the  altitude  of  the  aircraft  will  usually  offset  some  of 
these  disadvantages.  Dipoles,  slotted  monopoles,  and  surface 
wave  antennas  Cl4i  p. 29-30]  are  antennas  that  can  be  used  in 
airborne  applications. 

The  DF  function  of  the  system  could  use  any  of  the 
mentioned  types  of  antennas  if  the  system  uses  amplitude  or 
phase  comparisons  between  two  collectors.  C14:  p. 29-31] 
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Loop  antonnas  ar«  a  common  amplituda  sensitive  DF 
antenna  o-ften  used  in  a  system  called  a  goniometer  [1<S:  p. 
2h-14].  Horizontal  dipoles  or  dipoles  combined  Mith 
re-f lectors  are  also  antenna  types  that  are  used  in  DF  systems, 
as  well  as  Adcock  arrays  Ci6:  p.2h-i8]. 

Generally,  it  is  desired  to  have  high  gain  in 
antennas  -for  an  ESM  system,  and  the  speci-fic  application  will 
dictate  the  type  and  size  of  antenna  that  can  be  used.  More 
often  than  not,  an  omni-directional  antenna  of  some  type  is 
used  for  search,  and  a  rotating,  fixed,  or  multi-port  highly 
directional  antenna  is  used  for  direction  finding, 
b.  Receivers 

Most  discussions  on  EM  wi 1 1  address  the  receiver 
technology  used  to  detect  the  signal  brought  in  by  the  antenna 
system.  Several  types  of  receivers  are  in  use  by  the  Air 
Force,  Navy,  and  Army,  with  the  mission  usually  dictating  the 
choice  of  technology. 

A  quick  and  general  overview  of  receiver  types  and 
technologies  follows  in  this  section.  Although  some  of  the 
technologies  have  been  around  for  some  time,  and  some  still 
yet  to  come,  they  are  increasingly  being  dominated  by  the 
processing  requirements  and  the  impact  of  microcomputers. 

Cost  saving  and  simpl if ication  of  the  microwave  design 
requirements  are  a  result  of  this  increasing  influence.  In 
addition,  combinations  of  receiver  technologies  are  being  used 
to  overcome  limitations  of  individual  receiver  concepts. 


The  -future  holds  new  technologies  -for  the  detection  o-f 
signals  and  increases  of  speed  and  bandwidth  of  single 
receiving  devicesc  Surface  Acoustic  Wave  devices  will 
certainly  play  an  important  rule,  as  well  as  acousto-optic 
techniques.  Increases  in  computational  speed  and  capacity  may 
aid  or  deter  some  technologies  by  eliminating  some  of  the  need 
for  higher  sensitivity  in  the  detection  scheme. 

Table  2.2  C 171  gives  some  typical  ESM  receiver 
specifications. 


Table  2.2.  Typical  Receiver  Specifications. 


Frequency  Range  0.5  to  18  GHz 

Signal  type  100  ns  pulse  to  CW 

Sensitivity  <-70  dBm 

Resolution  3  MHz  (pulse) 

1  MHz  <CW) 

Ampl i tude  accuracy  1  dB 

Bearing  accuracy  3*  RMS 

Pulse  width  resolution  100  ns 

Time  of  Arrival  resolution  50  ns 
Signal  density  10^  pps 

Intercept  probability  169'/. 


< 1)  Direct  Detection.  Simplicity  in  design  gives 
the  direct  detection  receiver  advantages  in  cost,  reliability, 
size,  and  weight.  This  design  is  also  proven,  currently  in  use 
in  the  front  end  of  radar  warning  receivers  (RWR) ,  giving  a 
high  Probability  of  Intercept  (POD  for  signals  above  its 
Minimum  Discernabte  Signal. 

Today's  high  density  environments  pose 
problems  for  the  DO  receivers,  degrading  their  performance 
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quickly.  Frequency  resolution  is  not  very  high  as  the 
bandwidth  is  usually  an  octave  or  more  in  width.  C18:  p.26] 
Sensitivity  is  a  problem,  although  RF  ampl  i-f  ication  can 
provide  improvements. 

Quite  susceptible  to  ECCM,  the  DO  receiver 
cannot  handle  PRI  agile  signals,  and  does  not  provide 
-frequency  measurement,  an  important  sorting  parameter. 

Figure  2.2  is  a  direct  detection  receiver. 


Th«  operation  of  the  discriminator  is  the 


heart  of  the  receiver.  Note  Figure  2.3<a) .  Input  signals 
are  divided  into  two  parts  and  phase  correlated  after  one 
signal  is  delayed  by  a  known  time.  The  length  of  the  delay 
line  is  related  to  the  amount  of  frequency  ref>o1ution  desired 

I 

and  the  wavelength  that  corresponds  to  the  desired  bandwidth. 
By  banking  several  discriminators  and  different  lengths  of 
delay  line,  as  shown  in  Figure  2.3<b) ,  a  desired  frequency 
band  can  be  covered.  This  approach  allows  the  use  of 
practical  and  relatively  inexpensive  delays  rather  than 
attempting  to  use  one  delay  line  precise  enough  to  handle  the 
desired  frequency  region.  The  length  of  the  longest  delay 

becomes  the  frequency  resolution  of  the  bank.  CIS:  pp. 32-331 

( 

Measurement  of  the  frequency  of  the  incoming 
signal  with  these  analog  discriminators  <or  current  technology 
Digital  IFM^s)  allows  this  important  parameter  to  simplify  the 


sorting  process, 


On  the  down  side,  the  IFM  doesn^f  handle  high 


density  environments  very  well  unless  a  subsystem  that 
determines  when  2  pulses  occur  simultaneously  is  included. 
Also,  the  delay  lines  require  special  attention  <constant 
temperature)  and  this  usually  increases  the  required  input 
power  to  the  system. 

Morking  in  conjunction  with  an  analysis 
receiver,  the  IFM  makes  an  ideal  acouisition  receiver  by 


provi.ding  th«  -frequency  o-f  s.  signal  o-f  interest.  Analysis 
receivers  of  another  type  can  then  be  tuned  to  the  frequency, 
C19;  p. 32-345  . 

..  '  /  ,  •  '  i  ‘ 
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(3)  Superheterodyne.  Most  comnon  of  the 
receivers  in ; current  use*  the  superheterodyne  provides, hxgh 
s41ectivityi  good  rest stanco  to  Jamming ,  and  is  a  proven 
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DraufbacHs  include  sVom  time  in  scanning  the 
frequency  band  q*  interest  and  inability  to  see  frequency 

i  ‘  '  M  r 

agile  signals.  C19;  p.4S93  ■  \  - 
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Figure  2.6.  Superheterodyne  Receiver 


<4)  Channel ized.  Efforts  to  increase  the  POI  of 
superheterodyne  receivers  and  deal  v^ith  an  inability  to  handle 
high  density  environments  led  to  the  concept  of  a  bank  of 
filters  followed  by  individual  signal  detectors  that 
determine  when  a  signal  is  within  the  filter  bandwidth.  The 
RF  spectrum  is  simply  broken  into  pieces.  See  Figure  2.5. 


A  key  advantage  o-f  the  ihannelized  receiver  is 
removal  o-f  pulse  overlap  o-f  near  siniu  1  taneous  signals  (only 
pulse  -frequencies  that  exceed  the  channel  spacing).  C17: 
p.  les] 

By  -following  an  initial  bank  o-f  filters  with 
multiplexers  and  -fixed  oscillators,  then  repeating  the 
scheme  with  a  smaller  bandwidth  -filter,  the  RF  input  can  be 
down  converted  to  a  -final  baseband  -frequency  range  with  a 
higher  -frequency  resolution. 


Figure  2.S.  Channelized  Receiver 

<5)  Compressive.  This  receiver  can  be 
characterized  as  a  -fast  scanning  superheterodyne  receiver. 
Ideally,  the  local  oscillator  scans  the  RF  bandwidth  being 
covered  in  a  time  less  than  the  narrowest  pulse  to  be 
intercepted.  Its  high  POl  and  the  ability  to  handle  wideband 
signals  and  -frequency  agile  signals  makes  this  receiver  a  good 
choice  -for  an  ESM  receiving  system.  It  has  an  excellent 
ability  to  separate  signals  closely  spaced  in  -frequency. 
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Basically,  ths  compressiv*  rscaivsr  provides 
an  IF  signal  that  is  up-chirped  in  frequency  from  the  RF 
input.  The  compressive  filter  operates  in  such  a  May  that  the 
low  frequency  components  of  the  IF  signal  are  delayed  longer 
than  the  high  frequency  components,  with  the  result  that  the 
shape  of  the  output  pulse  of  the  detector  is  in  essence  the 
Fourier  Transform  of  the  IF  signal  pulse.  The  maximum 
amplitude  of  the  detector  output  falls  at  a  point  in  the  scan 
time  that  is  proportional  to  the  frequency  of  thw  RF  input. 
tl8:  pp.  38-311 

Criticality  of  the  alignment  between  sweep  and 
compressive  delay  makes  this  complex  technique  difficult  to 
manage  C19t  p.  4591.  But  it  is  at  the  forefront  of  some  of 
the  newest  technologies,  with  S^'s  performing  the  delay 
necessary  for  the  up  and  down  chirp.  Note  Figure  2.6. 


Figure  2.6.  Compressive  Receiver. 

<6)  BraoQ  Cel  1 .  An  acousto-optic  Bragg  cell  will 
interact  a  propagating  acoustic  wave  with  an  optical  beam  to 
give  a  diffracted  output  proportional  to  the  frequencies  and 


pow*r  pr«««nt  in  th*  RF  input  signal .  Th«  acoustic  wave, 
which  varies  the  index  o-f  re-fracticn  of  the  cell  material,  is 
generated  by  applying  the  RF  input  to  an  acoustic  transducer 
in  the  Bragg  cell.  A  laser  beam  is  deflected  and  modulated  in 
intensity  by  the  index  variations  which  are  proportional  to 
the  RF  input  frequency.  A  Fourier  Transform  lens  collects 
the  proportionately  diffracted  light  which  then  falls  on  a 
photodetector  array.  The  output  of  the  detector  array  is 
similarly  proportional  to  the  RF  input  f reqire.-tcy .  C  18i  p.28] 
Having  pluses  for  high  sensitivity,  high  POI ,  and 
high  selectivity  makes  this  receiver  a  choice  of  the  future. 

It  doesn-'t  handle  frequency  agile  sources  well,  however, 
because  of  the  time  needed  in  determining  the  output  of  the 
photodetector  array.  Pulse  width  measurement  is  also 
eliminated  with  the  use  of  the  detector  array  because  the 
array  only  responds  to  the  energy  contained  in  several  pulses 
from  the  lens.  C 19:  p.460]  Figure  2.7  is  the  receiver. 


Figure  2.7.  The  Bragg  Cell  Receiver. 
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2.  The  Ppoc«ssinQ  System 


A  steady  stream  o-f  pulse  descriptor  words  from  the 
receiving  system  must  be  analyzed  -for  the  vital  in-formation 
that  is  contained  in  them.  The  processing  system  must  extract 
this  in-formation  and  maKe  decisions  regarding  the  emitters 
whose  signals  it  has  processed.  In-formation  extraction, 
analysis  and  identi-f ication  must  be  -fast  to  avoid  a  mission 
and  li-fe  threatening  situation. 

This  section  will  discuss  the  overall  purpose  and 
general  operation  o-f  preprocessors  and  processors  in  the  ESM 
system.  A  more  complete  discussion  will  be  -found  in  Re-ference 
28. 

a.  Preprocessor 

Basically,  the  preprocessor  must  prepare  the 
signal  received  from  the  receiving  system  for  preliminary 
analysis  and  for  advanced  analysis  in  the  main  processor. 
Figure  2.8  is  an  ESM  preprocessor.  C20:  p.81 

After  digitization,  the  initial  parameters  of 
pulse  width,  time  of  arrival,  direction  of  arrival,  amplitude, 
and  frequency  must  be  measured  on  a  pulse  by  pulse  basis. 

(See  Parameter  Measurements  in  this  chapter.) 

The  pulse  data  must  be  compared  with  previously 
received  and  analyzed  signals  to  see  if  this  signal  has 
already  been  processed.  Pulse  data  is  compared  with  signal 
data  in  the  files  that  contain  the  active  emitters  and  the 
uninteresting  emitters.  If  a  match  is  present,  then  the  data 
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can  b«  d»l«t«d  -froni  th«  procass*  Reduction  o-f  th«  data  atraam  j 

3 

! 

and  sorting  o-f  th*  signals  is  tho  purposo  and  rosul  t  of  this 

dalstion  in  addition  to  avoiding  unnocossary  procassing.  Ona 

of  tha  main  Jobs  of  tha  praprocassors  is  to  convart  tha 

approx imatal y  10 ■  to  10*  pulsas  par  sacond  baing  racaivad  to 

about  1000  pulsas  par  sacond  that  a  good  main  procassor  can  ] 

handla.  C21t  p.l<S4] 

Unmatchad  or  naw  signals  must  ba  daintarl aavad 
Xsaa  daintarl aaving  in  this  chaptar)  and  procassad  in  tha  PRI 
procassor  saction.  Calculation  of  the  PRI  is  usually  dona 
simply  by  subtracting  TOA'^s  of  similar  pulsas.  Staggarad 
PRI's  and  Jittarad  PRIG'S  have  to  ba  handled  with  a  more 
complex  algorithm.  C20:  p.  121 

Mith  a  Known  PRI,  TOA  prediction,  additional 

J 

sorting,  and  comparison  of  calculated  PRI's  is  dona  until  a 
PRI  is  confirmed  to  exist  in  tha  data.  This  information,  tha 
newly  calculated  PRI  and  tha  initial  parameters  of  tha  pulse, 
is  than  passed  to  main  memory  of  tha  procassor  for  analysis 
and  updating  of  tha  active  files  so  that  pulsas  of  tha 

1 

incoming  stream  can  ba  delated  as  soon  as  possible. 

i 

In  summary,  tha  preprocessor  prepares  tha  incoming  j 

signals  for  processing,  measures  their  initial  parameters, 
performs  preliminary  analysis  and  sorting,  and  passes  tha  data 
information  to  tha  main  procassor. 
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b.  Main  Processor  C20s  pp. 13-183 

After  all  the  action  in  the  preprocessor  and  the 
data  reduction  that  has  been  accomplished,  there  is  still 
plenty  to  do.  The  functions  of  the  main  processor  arex 

1.  Advanced  analysis  of  the  signal  for  scan  rate  and  exotic 
emitter  (such  as  chirp  radars)  identification 

2*  Emitter  and  platform  identification 

3.  Updating  the  emitter  files 

4.  System  control  and  monitoring 

9.  Interfacing  uxith  the  operator  and  ECM  equipment 

Since  the  processor  is  usually  a  general  purpose 
computer,  additional  functions  of  navigation  or  routine  system 
calculations  could  be  added.  It  can  also  help  the  PRl 
processor  handle  its  calculations. 

Emitter  identification  is  accomplished  by 
comparing  the  measured  emitter  parameters  with  emitter 
parameter  lists  until  an  identification  is  made  or  the  signal 
is  tagged  as  unidentified.  Interface  with  the  operator  for 
manual  analysis  is  often  the  result  of  an  unidentified 
signal . 

3 .  Display  System 

An  operator  must  have  a  visual  presentation  of  the 
signal  and  its  parameters.  Almost  universally  CRT  displays 
are  used.  CRT's  usually  have  an  intensity  modulated  or  a 
deflection  modulated  design,  with  advantages  and  disadvantages 
in  each  one. 

Multi-gun  CRT  displays  can  create  a  combined  display 
that  offers  a  pulse  analyzer,  presenting  frequency,  PRF,  PW, 
and  the  presentation  of  an  amplitude  spectrum. 


A  wjLrning  ftystem  u«uai11y  op«rAt*«  in  conjuction  with 
detection  o-f  high  threat  signals.  Flashing  lights  and/or  an 
aural  alarm  will  occur,  indicating  immediate  attention  is 
required. 

The  audio  system  operates  in  parallel  with  the  warning 
system.  0-f-fering  not  only  alarms  to  high  threat  signal 
reception,  the  system  usually  has  accomodations  that  allow  the 
operator  to  monitor  the  received  signal  through  headsets  or  a 
speaker.  A  trained  and  experienced  operator  can  determine 
scan  rate  and  analyze  the  scan  modulation  type  i-f  the  signal 
density  is  low  enough. 

Bearings  -from  DF  equipment  are  usually  presented  on  a 
polar  display.  The  direction  is  usually  indicated  in  relative 
bearings  oi  0  to  360*  around  the  plat-form,  with  some  output  or 
means  o-f  determining  the  true  bearing  -for  plotting  purposes. 

Computer  management  o-f  an  CSM  system  allows  -for 
most  calculations  o-f  bearing,  FRF,  PM,  and  frequency  to  be 
accomplished  and  shown  as  soon  as  the  bearing  is  displayed  to 
the  operator.  An  excellent  example  of  computer  managed 
displays  is  the  ALQ-'78  ESti  system  on  the  U.S.  Navy  P3-C 
Orion.  Digital  readouts  of  the  PRF,  PM,  and  frequency  are 
available  next  to  the  displayed  bearing  from  the  aircraft. 

The  bearing  remains  on  the  screen  and  updates  automatically 
until  removed  by  the  operator  or  the  system.  Mith  several 
cross  bearings,  the  emitter  location  can  be  obtained  since  the 
computer  is  not  only  managing  the  display  but  the  navigation 
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of  th«  Aircr&ft.  By  "hooking"  the  place  where  the  bearings 
cross,  the  operator  can  cause  the  computer  to  calculate  the 
1  at/long  of  the  emitter  fix. 

B.  ESM  SIGNAL  PROCESSING 

The  acquisition  of  signals  is  usually  a  design  problem 
that  is  driven  by  the  characteristics  of  the  signals  to  be 
acquired.  The  designer  chooses  antennas  and  receivers  to 
detect  signals  with  regard  to  the  mission,  the  platform,  and 
the  frequency  band  of  interest.  Generally,  he  desires  the  ESM 
system  should  have  rapid  acquisition,  high  sensitivity  and 
frequency  resolution,  and  wide  frequency  coverage. 

After  these  signals  are  acquired,  signal  processing  is  the 
order  of  the  day.  Processing  is  required  immediately  to 
determine  the  high-threat  signals  <especially  for  tactical 
aircraft) ,  and  a  lot  of  processing  is  required  in  todays s 
dense  electromagnetic  environments.  The  signals  are  getting  ‘ 
more  complex  also,  with  Jittered  and  staggered  PRF^s  and 
variable  PM'^s  in  addition  to  frequency  agility  and  new 
schemes  of  radiating  and  detecting  a  return. 

Computational  speed  becomes  one  answer  to  solve  these 
processing  riddles,  and  the  debate  between  digital  and  analog 
techniques  fuels  upon  this  quest  for  speed. 

Surface  Acoustic  Mave  (SAM)  devices  and  Charge  Coupled 
devices  (CCD's)  are  being  challenged  by  digital  technology 
staked  by  the  Pentagon^'s  Very  High  Speed  Integrated  Circuits 
Program  (VHSIC) .  These  digital  IC's  are  being  structured  to 


perform  high  spMd  ca.1  cul aition ,  with  a  N  X  N  matrix  inversion 
<a  operation)  being  a  gauge  o-f  successful  application. 
Today's  analog  technologies  currently  offer  high  speed  also, 
but  will  require  more  precision  and  dynamic  range  to  meet  this 
challenge.  C22:  p.911 

In  the  end,  the  ESM  system  designer  will  win  as  new 
d\gita1  architectures  and  new  device  structures  will  yield 
increase  in  speed  and  throughput. 

This  section  will  discuss  important  signal  processing 
parameters  of  the  ESM  system,  and  the  techniques  used  to 
separate  the  signals  for  identification  and  processing  from 
the  imnense  amounts  of  pulses  being  received.  The  entire 
process  is  a  continuous  flow  of  calculations,  comparisons, 
and  decisions  involving  tolerances  in  frequencies,  pulse 
width's  <PW> ,  pulse  repetition  intervals  (PRl) ,  and  other 
parameters. 

1 .  Parameter  Measurement 

Listed  below  are  the  initial  parameters,  as  provided 
by  the  receiving  system  to  the  preprocessor .  They  are  the 
basic  tools  that  are  used  to  dein ter leave  <see  next  section) 
the  pulses  and  match  them  to  existing  tracked  or  not  tracked 
emi tters. 

1.  Direction  of  Arrival 

2.  Amplitude 

3.  Frequency 

4.  Pulse  Midth 

5.  Time  of  Arrival 
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Those  thst  sren't  matched  must  be  analyzed  further  by 
ca1culatin9  a  pulse  PRI  and/or  undergoing  advanced  analysis. 
Calculation  of  the  PRI  is  done  in  the  preprocessor. 

Advanced  analysis  is  accomplished  in  ^he  main  processor  using 
derived  higher  order  parameters  during  the  tasK. 

A  discussion  o-f  some  aspects  of  how  the  initial 
parameters  and  higher  order  parameters  are  measM'*®d  follows  in 
this  section.  121:  p.l66] 

a.  Direction  of  Arrival  <00A> 

An  important  parameter,  the  DOA  could  be  measured 
with  a  rotating  OF  antenna  but  this  greatly  lowers  the 
Probability  of  Intercept  <POI> .  A  more  usual  method  uses  a 
multi-port  <antennas>  bearing  measurement  subsystem  that  has  a 
high  POI  . 

Phase  comparison  techniques  can  also  be  used, 
providing  greater  accuracy  but  with  a  tradeoff  of  smaller 
bandwidth . 

b.  Amplitude  ^ 

Usually  the  amplitude  is  termed  by  the  peak  value 
of  the  received  pulse  with  some  tolerances  included  to  provide 
for  variances  in  emitter  radiated  power. 

The  amplitude  parameter  is  used  in  the  ad'.'anced 
analysis  section  of  the  typical  processor  to  determine  a  scan 
pattern,  but  is  not  considered  a  good  sorting  parameter 
because  of  propagation  corruption,  multi -path,  and  other 
variables  that  influence  the  amplitude  of  the  received  pulse. 
C16:  p.  2g-3I 
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■'  c.  Fr«ciu«ncy  '  '  i  -  .  >  .  , 

Fr«qu*ncy  may  be  measured  with  a  Ssianning  < 

‘  ‘  '  i  ...  f  i 

superheterodyne  receiver,  but.  this  method  has  a  <  ow  POl  -for 

;  ■/  •  ,  ■  ■ ,  • ,  ,■ 

reasons  similar  to  those  of  using  a  rotating  antenna  -for  DF. 

A  channel  ized  receiver  provides  a  f requ measurement  but 

■('  ' 

this  usually  4s  a  much  more  expensive  system.  ■  i 

'  <■-  ‘  ■- 'J  '  - 

Most  common  is  the  IFM  receiver  <see  section  on 
receivers  in  this  chapter).,  which  has  an  excellent  resolution 

'  -  ■  ■  .  -  ■  -  -  y,  V  \ 

^  and  high  POI  .  .I 

d.  Pulse  Width 

'  ■  I  .  i  \ 

Pulse  widths  can  be  measured  by ^simply  noting 
when  the  pulse  rises  and  falls  through  a  threshold. 

Reflections  tend  to  severely  corrupt  this  method  o-f  ' 
measurement,  howerver . 

^  { 

A  more  accurate  method  is  to  sample  the  pulse  many 

times  and  take  the  PW  in  the  distance  between  the  3  dB 
points.  , 

,e.  Time  of  Arrival  <TOA) 

1 

Threshold  crossing  measurement  of  this  parameter 
could  also  be  used  as  in  the  case  of  the  PW,  but  a  better 
measurement  can  be  obtained  by  noting  the  TOA  of  the  first  3 
dB  point.  TOA^s  are  used  in  the  calculation  of  PRI^s. 

Adjacent  pulse  TOA^s  are  simply  subtracted  to  find  the  PRI . 
f .  Higher  Order  Parameters 

Parameters  that  are  calculated  from  the  other 


initial  parameters  are  the  derived  or  higher  order  parameters. 
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The  PRI  is  used  to  link  interrupted  or  split  . 
chains  of  pulses  from  the  same  emitter.  It  can  also  be  used 
as  a  parameter  to^  identify  a  Jittered  or  staggered  PRI . 

Scan  patterns  can  be  recognized  by  analysis  of  the 
pulse  amplitudes  generated  as  the  radar  beam  sweeps  across  the 
platform.  Knowledge  qf  the  scan  pattern  can  be  cf  help  in 
id^fotifying  the  emitter  since  the  radiatiqn  pattern  can  be 
correlated  to  emitter  type.  C21:  p. 1671 

^  /'’  I  ■  '  i.  '  .  .  j  - 

I  .  /  '  '  *  - 

2.  Oeinterl  eaving 

A  stream  of  pulses  from  several  different  emitters 

'  f  ^  ■ 

^  includes  different  frequencies,  amplitudes,  and  PRI-'s  mixed 
together  < inter  1 eavedi .  Separating  these  different  pulses 
into  chains  of  similar  pulses  for  identification  is  the 

■  ■  ■  ■  '  ,  f  ,  -  ' 

deinterleaving  problem. 

With  high  density  environments,  the  time  available  for 
processing  each  pulse  may  be  only  around  1  microsecond.  C21> 
p.l68]  Therefore,  deinterleaving  algorithms,  if  done 
digitally,  must  be  fast.  Analog  equipment  that  sort  pulses 
have  the  same  requirement.  This  required  speed  is  the 
obstacle  that  hinders  the  use  of  Digital  FFT‘'s  for  pulse 
sorting. 

Provision  must  be  made  in  algorithms  for  missing 
pulses  in  the  data  stream.  Missing  pulses  is  one  of  the 
greatest  confusion  factors  in  the  system.  Methods  are 
currently  in  use  to  smooth  fluctuations  of  the  number  of 
pulses  that  are  input  into  the  preprocessor.  C21i  p.l671 
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General  dis-cussions  of  deinterleaving  methods  are 
given  in  this  section.  The  first  two  are  digital  methods, 
using  the  measured  parameter  of  the  pulse,  and  use  of  one  or 
the  other  depends  upon  the  accuracy  of  the  ESH  system.  Two 
hardware  methods  are  briefly  covered  also. 

a.  Ceil  or  "Pigeoti  Hole"  Techniques 

A  systm  that  accurately  measures  the  parameters 
of  the  pulse  stream  can  use  this  technique  as  a  fast  and 
efficient  deinterleaving  scheme. 

Pulse  data  with  similar  parameters  are  directed 
into  cells  or  "pigeon  holes*.  The  cells  soon  contain  alike 
parameters,  such  as  bearing,  frequency,  or  pulse  widths  that 
can  be  analyzed  for  emitter  identification  or  further  sorting 
if  too  wide  of  a  variation  exists. 

Variations  could  be  caused  by  emitter  variations 
itself  or  purposeful  variations  such  as  a  jittered,  staggered 
PRI ,  or  frequency  agile  radar.  If  variations  (unmatched 
parameters)  are  present,  then  they  are  passed  to  a  more 
complex  algorithm  for  processing  for  the  special  modulations. 
Otherwise,  a  simple  analysis  or  comparison  with  the  active 
emitter  file  or  emitter  parameter  file  results  in  an  emitter 
identification . 

b.  Time  of  Arrival  Techniques  C21:  p.l68] 

When  accurate  computation  of  pulse  parameters 
isn^t  available,  a  time  slice  of  the  pulse  stream  is  used  to 
simplify  the  data  under  analysis.  A  time  slice  will 


contain  about  8  to  512  pulsva.  This  schents  reduces  the  number 
o-f  pulses  being  considered  in  a  high  density  enviroment,  and 
forces  comparison  of  parameters  to  those  only  in  the  present 
time  si  ice. 

The  first  pulse  in  the  time  slice  is  used  as  a 
reference  and  the  rest  of  the  pulses  is  scanned  for  a  match. 
Mhen  a  match  is  founds  the  TOA  of  the  pulses  is  used  to 
calculate  a  PRI .  Mhen  a  sufficient  rumber  (usually  6  or  7)  of 
matching  pulses  are  founds  then  the  data  is  jrouped  as  a  chain 
and  is  output  to  the  main  processor.  A  problem  generated  with 
this  technique  is  the  splitting  of  chains  in  the  pulse  data 
stream  by  the  time  slice.  Several  time  slices  may  have  to  be 
analyzed  before  enough  matches  in  a  slice  are  found  to  form  a 
chain . 

c.  Analog  Methods 

A  pulse  sorter  operating  on  the  pulse  period 
th  has  some  analogy  to  digital  TOA  methods  is  currently  in 
use.  C16i  2g-'7  to  2g-?]  Consider  Figure  2.19. 

The  pulse  stream  is  fed  into  a  coincidence 
detector-discriminator  that  attempts  to  match  the  pulses  with 
gates  generated  by  a  gati.^g  pulse  generator.  The  gating 
pu'ses  are  initially  generated  at  a  PRI  slightly  higher  than 
the  highest  pulse  PRI.  By  gradually  reducing  the  gate  PRI, 
the  highest  PRI  pulse  train  will  be  captured  by  the  gates  if 
the  phase  difference  betweer  the  gates  and  pulses  of  the 
stream.  This  phase  difference  is  used  to  generate  an  error 
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signal  that  causas  tha  gata  PRt  to  dacraasa  until  synchronized 
with  the  highest  PRI  pulse  train  o-f  the  stream.  The  pulse 
train  is  now  sorted  out  o-f  the  stream,  and  the  remaining 
pulses  are  applied  to  a  second  coincidence 
detector-discriminator  where  the  process  is  repeated. 

Other  methods  are  used  on  the  pulses  left  in 
the  stream  after  passing  throught  all  the 

detector-discriminators.  Staggered  PRI's  can  also  be  handled 
by  using  two  or  more  gate  generators  operating  in 
synchronization . 

d.  Frequency  Domain  Methods 

Bandpass  comb  filters  Cl^:  p.2g-91  can  be 
used  to  pass  only  harmonically  related  frequency  components. 

By  tuning  the  filter  with  the  pulse  stream  as  an  input,  the 
filter  will  pass  components  when  it  is  tuned  to  a  PRI  of  the 
pulse  stream.  Disadvantages  of  this  method  include  the  low 
resolution  of  the  sorted  PRI's  and  inability  to  distinguish 
harmonic  components  of  different  PRI's. 

Digital  FFT‘'s  are  attractive  as  a  sorting 
tool ,  but  hardware  (size  and  cost)  and  computational  speed 
have  a  way  to  go  before  it  offers  practical  alternatives. 
V)H81C  will  certainly  improve  both  the  hardware  and  speed 
aspects. 
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III.  THEORY  OF  kMLSH  AND  RAOEMACHER  FUNCTIONS 


Electrical  Engineering  could  be  said  to  revolve  around  the 
sine  and  cosine  -functions.  They  are  the  basis  -for  development 
in  many  areas  due  in  part  to  the  inherent  properties  in  their 
-frequency  domain  representations. 

The  application  of  digital  techniques  and  semiconductor 
technology!  though i  has  brought  -forth  uses  and  awareness  of 
other  orthogonal  functions.  These  often  do  not  have  the 
desirable  properties  of  sine  and  cosine  functions  for  use  in 
linear)  time  invariant  networks)  but  they  do  have  other 
advantages  that  render  them  useful  in  other  applications. 

The  early  part  of  this  century  saw  the  introduction  of 
several  two  level  orthogonal  functions,  or  sometimes  called 
binary  functions  because  of  their  amplitudes  taking  only  two 
values.  Work  in  orthogonal  matrices,  by  Sylvester  in  1867 
and  Hadamard  in  1893,  was  an  early  approach  to  these 
functions. 

In  1910)  Alfred  Haar  presented  a  complete  set  of 
rectangular  functions  that  took  only  two  values  but  yet 
provided  complete  expansion  of  a  continuous  function.  These 
Haar  functions  could  be  made  to  converge  uniformly  and 
rapidly.  See  Reference  23  for  additional  information. 

A  German  mathematician,  H.  Rademacher ,  presented  another 
set  of  orthogonal  functions  in  1922.  These  were  followed  by 
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I  -m  t  .  M  f  M  flLa  f. 


th«  W«i1sh  Functions,  d*-fin*d  in  1923  by  th«  Amsricsn 
ffluthsmsticisn  J.L.  Wsish.  Ths  RatdwRachsr  functions,  although 
indopondsnti y  prssontod,  w«r*  -found  to  b*  an  incomploto  but 
tru«  subsot  to  tho  Malsh  functions.  [23:  p.v-vil 

This  chaptor  will  briofly  oxamino  tho  Radomachor 
functions.  Thoy  aro  important  bocauso  products  of  th# 
functions  yiold  a  contain  ordoring  of  tho  Walsh  functions. 

Tho  Walsh  functions  will  thon  bo  oxaminod  in  mono  dotail. 

Tho  roador  is  roforrod  to  Appondix  B  for  a  gonoric  troatmont 
of  orthogonality  and  orthogonal  functions. 

A.  RADEmCHER  FUNCTIGNS 

Thoso  two^lovol  orthogonal  functions  aro  roprosontod  by 

R<n,t)  <3.l) 

n  ^  0, 1,2,3,... 

and  can  bo  soon  in  figuro  3.1  for  n  »  8  to  6.  C23i  p.6] 

Notico  that  thoy  aro  a  sorios  of  roctangular  pulsos  or 
squaro  wavos  and  havo  2*^^^  poriods  ovor  a  timo  baso  from  8 
to  T.  Mathomatical ly  thoy  could  bo  dofinod  by 

R<n,t)  -  sign(sin<2''lCt))  <3.2> 

Tho  first  Radomachor  function,  R(8,t>,  is  oqual  to  ono  for 
tho  ontiro  intorval ,  and  subsoquont  Radomachor  functions  havo 
odd  symmotry  with  amplitudos  of  -•’1  and  -1.  All  havo  unit 
mark-spaco  ratio. 

Ono  can  gonorato  thoso  functions  with  a  sino  function  of 
appropriato  froquoncy  with  amplification  and  hard  limiting. 
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Th«  appropriat*  -fraqucncy  moulci  b*  on*  with  th*  samo  z«ro 
crosaing  poaitiona  aa  tha  Radamachar  functions.  C23:  p.7] 


B.  VMLSH  FUNCTIGNS 

Of  mor*  importanca  than  tha  othar  two-taval  orthogonal 
functions  diacuaaad  ara  tha  Walsh  functions.  Thay  form  a  sat 
of  ractangular  wavaforma  taking  only  two  amplitudaa,  and 
-if  but  do  not  hava  a  unit  mark-apac*  ratio  lika  tha 
Radamachar  functions.  Thay  ara  dafinad  ovar  a  tima  intarval 
T,  which  must  ba  known  to  assign  valuas  to  tha  functions. 

Two  argumants  ara  raquirad  for  complata  spacif ication  of 
a  Wal sh  function 


MAL<n,t> 


<3.3> 


Th»  tim*  b««*  Ct]  o-f  th*  -function  is  usually  specified  as 
t/T,  and  thus  normalized  from  9  to  1.  The  number  Cnl  is 
equal  to  the  number  of  zero  crossings  a  Malsh  function  has 
during  the  time  base.  [23i  p.7] 

Figure  3.2  shows  the  Malsh  functions.  The  symmetry  of 
the  functions  and  the  concept  of  sequency  <next  section)  was 
used  by  Harmuth  (243  to  define  another  notation  for  each 
Ma1sh<n|t).  Each  MAL<n,t>  is  either  odd  or  even  about  the 
midpoint.  If  the  function  is  odd,  then  it  can  be  referred  to 
as  SAL(k,t).  If  the  function  is  even,  then  HAL(n,t>  is  also  a 
CAL<k,t) .  Note  the  similarities  between  this  notation  and  the 
sine-cosine  functions.  Relationships  between  Cn3  and  Ckl  are 

CAL<k,t>  «  li«^L<2k,t>  <3. 4) 

SAL<k,t)  »  WAL<2k-l,t> 

This  notation,  SAL(k,t)  and  CAU<k,t>,  is  seen  on  the  right 
side  of  Figure  3.2. 

i.  Sequency,  Ordering,  and  Phasing 

Sequency  is  a  term  proposed  by  Harmuth  to  describe  a 
periodic  repetition  rate  which  is  independent  of  waveform.  It 
is  defined  as  "one  half  of  the  average  number  of  zero 
crossings  per  unit  time  interval*  and  abbreviated  *Zps* .  C23: 
p.l3]  A  109  HZ  sine  wave  has  209  zero  crossings  per  second, 
so  frequency  is  a  special  measure  of  sequency  as  applied  to  a 
sinusoidal  waveform. 

Ordering  the  functions  by  ascending  number  of  zero 
crossings  (MAL<n,t>  notation)  is  called  sequency  ordering. 
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WA».IZ9,T|_  _  _  _  _  __  „  »  M  „  „  n  <— CaL(i4,T) 


WAL(28,T)“ 
WAL(27.T)  - 

W*L1Z6,T)  - 
WAL(2S,T)- 
WAL(24,T)” 
WAl.(23,Tl- 

WAL(22,T)- 

WAL(ZI.T)- 

WA1.(20,T)“ 

WAU(I9.T»- 

WALOe.l )  - 
WAC  (I7,T). 
WAL<I*,T)- 

WAtdS.T)- 

WAt(l4,T)- 
V>AL(I9,T). 
WAL  (12, T)  ' 

WAU  (ll,T)- 

WALdO.T). 

WAL  (9.T)  . 
WAL  (B.T)- 
WAL  (7,T)- 

WAL  (•,T)' 
WAL  (9,T)  , 
WAL  (4,T). 

WAL  (3,T). 


WAL(2,T)- 
WAL(I,T). 
WAL  (0,T)- 


-SAL(!4,T) 

-CAL(I3.T) 
■SAL1I3,T) 
-CAL  (12, T) 

-SAL(I2.T) 
.CAL  (II, T) 
-SALdI  ,T) 
.CAL  (10, T) 

-SAL(lO,T) 

,CAL(9,T* 

-8AL(9,T) 

-CAL(8,T) 

-SAL  (8.T) 
-CAL(7,T) 
"SAL(7,T) 
-CAL  (8,T) 

~SAL(e,T) 

.CALlS.T) 

-SAlI5,T) 

_CAL(4,T) 

-  SAL  (4,T) 

-  CAL  13, T) 

-  sal  (3,T) 

-  CAL  12, T) 


•SAL(2,n 
.  CALd.T) 

■sal(i,t) 

-CAL(0,n 


Figur*  3,2.  The  Walsh  Functions  in  Sequency  Order 


-'a  r',’  ^  /Vi- A '  -.C 


Although  Cn3  is  not  the  sequency  of  the  function,  the 
equivalence  o-f  the  WAL<n,t)  uiith  3AL<K,t>  and  CAL<K,t> 
notations  automat' cal  ly  arranges  the  -functions  in  ascending 
order  o-f  sequency.  The  value  o-f  Ck]  in  the  SAL  and  CAL 
-function  is  the  sequency  -for  the  -function.  Seq^^ency  ordering 
is  the  pre-f erred  ordering  -for  spectral  analysis  and  -filtering 

The  natural  ordering  <or  Pal ey  ordering)  is  obtained 
by  generating  the  Walsh  functions  with  products  of  Rademacher 
functions.  Figure  3.3  shows  the  Malsh  functions  in  natural 
order,  which  are  referred  to  as 

PAL€n,t)  <3.S) 

n  ^  0il,2f3|... 

because  Paley  C24]  used  this  ordering.  For  theoretical  and 
mathematical  work,  image  transmission,  and  computational 
efficiency,  this  ordering  is  usually  the  preferred  one.  C23: 

p.  181 

A  modulo**2  addition  <Appendix  C)  relationship  exists 
between  the  Walsh  functions  with  sequency  ordering  and  the 
functions  with  natural  ordering.  C23i  pp. 31-321  An  example 
Mill  illustrate  the  conversion  between  the  two  orderings. 
Consider  WAL<n^,t)  «  WAL<13,t). 

Here  n^,,  »  13  *  110l2<  To  find  the  equivalent 
natural  order  function,  PAL<np,t>,  first  add  a  zero  <base 
2)  as  the  leftmost  digit  to  base  2  n^.  Then  generate  the 
value  of  np  in  base  2  by  doing  modulo-2  addition  between 
consecutive  zeroes  and  ones  of  base  2  n^,  starting  with  the 
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.'I'.'T.'.'iv'-. 


P*L  (31,1) . 
PAUOO,!) 


•  WALCa.t) 
■  V«IL(ZO,t) 


PM.(29kO- 
PAL  (28,1)  - 


■  WA(.(22,I) 
-  WAU23,>) 


PAI.(2T,t). 
PAL  (28.1)  • 


.  WAL(l8,t) 


PAL(23,I)- 
PAL  (24,0' 


.  WAL(I9,() 
-  WAL(I7,|) 


■  WAL(I6,<) 


PAL  (23,0- 
PAL(22,0  ■ 


'  WALUS.t) 


PAL  (21,1)  . 
PAL  (20,0- 


.  WAL{27,i) 
•  WAL(25,I) 


PAL  (19,0  • 
PAL  (18,0 ' 


■  WAL(24,t) 
.  WAL(29,0 


WAL(28,I) 


PAL  (.7,.)  WAL(S0.t) 

WAL(3I..? 


PAL  (13,1) 
PAL(I4,0 


PAL  (131). 
PAL  (t2,t)  - 


pal  (M,0- 

PAL  (10,0  ' 


WAL(3l,t? 

WAUIO.I) 

WAL(II,I) 

WAL(9,I) 

«AL(8,0 

WAL(I3,0 

WALdS.t) 


pal  (9,0- 
pal  (8, II 


.  WAL(I4,I) 


PAL  (7,0  , 
PAL  (8,0 


.  WAL(I3.I  ) 
•«Al(9.i) 


■  W»L(4,l) 


pal  (9,1). 

PAL(4,t|' 


PAL  (3,1)- 
PAL  (2,0  ■ 


■  WAL(8,I) 

-  WAl(7,I) 
.  WAL(2,t) 


PAL  (1,0- 
PAL  (0,0  • 


.  WAL  (3.t) 
■  WAL  (1,1) 

.  WAL  (0,0 


Figur*  3.3,  Th«  M«1 «h  Functions  in  Natural  Order. 
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l«»<ftmoftt  digit  <  the  zero  th«t  was  added  on  to  base  2  n^>  . 

Following  these  rules,  *  13  =*  11812  =®  0110l2»  and 

0  •  1  »  1  a  k3 
1  •  1  *  0  as  k2 

1  •  0  *  I  a  kj 
0  •  1  s  1  «  k0 

So, 

np  a*  ^k3k2kjK0)2  *  10 1 12  **11 

Figure  3.3  shows  the  result  is  correct. 

!r>WL<  I3,  t)  -  PAU(  ll,t) 

The  phasing  o-f  the  -functions  as  given  in  Figure  3.2 
and  3.3  is  known  as  Harmuth  C2S]  phasing.  This  phasing 
emphasizes  the  phase  similarities  with  sine-cosine  -functions. 

t 

Positive  phasing  C23:  p.l83  is  where  all  of  the 
functions  begin  at  -*-1.  This  involves  a  sign  change  for  many 
of  the  Figure  3,2  and.  3.3  functions,  and  is  the  result  of 
function  derivation  from  Hadaraard  matrices  (orthogonal 
matrices  composed  only  of  -*-1  and  -1  elements;  see  Ref.  23, 
pages  24-2S)  or-  Rademacher  functions  (next  section)  . 

2.  Derivation  of  the  Walsh  Functions 

Products  of  selected  Rademacher  functions  will  yield  a 
complete  set  of  Walsh  functions  in  natural  order.  Recall  that 
the  natural  ordering  of  the  Walsh  functions  was  referred  to  as 
P#^L<n,t)  . 

In  terms  of  Rademacher  functions, 

m■^  1 

PAL<n,t)  -'TX'  <3. 6) 

i»l 

where  is  either  a  zero  or  a  one,  indicating  the  presence 
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,or  not  the  i'th  Radecnacher  -function.  To  find  the  bi's 
and  Cm!S ,  represent  Cnl  as  a  binary  number, 

u#i  th  m  »  the  highest  poMer  o-f  two  found  in  the  binary  number 
jlnl.  The  b^'s  present  in  the  positions  of  base  2  Cnl 
indicate  the  presence  or  not  of  the  Rademacher  function  in 
Equation  <3,4> ■  ,  , 

Consider  PAt,<9,t)  .  Then 

\  l  ay' 

n  »  1001 2 

with  m  s  3,  tj^  =  1,  bj  =  1,  and  finally, 

PAL<9,t)  »  R<4,t>R<  l,t) 

\ 

It  should  be  noted  that  this  product  is  positive 
phasing.  Harmuth  phasing  can  be  obtained  by  defining  the 
Rademacher  functions  over  the  interval  This 

definition  inverts  all  of  the  functions  shown  in  Figure  3.1. 
Products  of  the  Inverted  functions  will  yield  natural  ordered 
and  Harmuth  phased  Malsh  functions.  123:  p.22]  The  Malsh 
functions  with  sequency  ordering  can  be  derived  with 
Rademacher  functions  by  use  of  the  C^ray  Code  (Appendix  C>  . 
Consider  MAL(13,t>  where  n  «  13. 

In  Gray  Code,  n  13  is  represented  by  1011,  and 
recognizing  the  ones  of  the  Gray  Code  are  in  the  fourth, 
second,  and  first  positions,  (similar  to  the  derivation 
of  the  PAL(n,t)  functions),  C23:  p.23] 

WAL(l3,t)  «  R(4,t)R(2,t)R(  1,  t> 


other  derivations^  oi  the  Walsh  -functions  can  be 


accomplished  by  dif-fsrence  equations,  Hadamard  matrices,  and 
doolean  synthesis.  See  Reference  23,  pages  20-26. 

3 .  Walsh  Series 

It  is  iMell  known  that  a  time  function  f(t>  can  be 

expressed  as  a  sum  of  a  series  of  sine  and  cosine  functions. 

/  * 

Each  function  has  a  coefficient^  that  determines  the  value  of 
the  function  in  that  series. 

The  time  function  can  be  expressed  in  a  similar  way 
using  the  Walsh  functions.  C23:  p.l3] 


f<t)  a  a0l«ML(0,t>  *  Z  a„WAL(n,t> 

n“l 

where  the  coefficients  would  be  calculated  from 


<3.8> 


a0/2  *  1/T  •  Cf<t>l«WL<0,t>dt 

J0 

an  -  l/T  .  ^f<t)WSiL<n,t)dt 


<3. 9) 


<3. 10) 


with  N  »  the  desired  number  of  Walsh  terms,  and  T  is  a  given 
length  of  the  function. 

Each  Walsh  function  has  a  coefficient  associated  with 
it  that  gives  the  value  of  the  function  in  the  series  the  same 
way  as  does  the  Fourier  Series  coefficient.  [23:  p.l3] 

If  the  function  is  periodic,  and  T  is  normalized  to  be 
equal  to  1,  then  a0/2  is  the  mean  of  the  function,  since 
WAL(0,t)  is  equal  to  1  over  the  interval  0it<l.  The 
expressions  for  a0/2  in  both  the  Walsh  and  Fourier  Series 
are  equal  with  this  normalization. 


Consider  that  the  -function  is  absolutely  integrable 
in  the  interval  8<t<l.  Then  the  -^unction  can  be  expanded  in  a 
Malsh  series  o-f  the  form  C26:  p.232  and  23:  p.40] 

m 

f<t>  *  Z  a(n>l»WL<n,t)  <3.  11) 

n«»0 

where 

a<n)  *  ^S<t)l»h^L<n,t>dt  <3.12) 

A  more  general  definition  for  a  periodic  function  defined  over 
an  interval  <0,T)  is 

w 

f<t)  =  r  a<n>l*WL<n,t/T)  <3.1U) 

nai0 


acn)  «  1/T  T  f<t)WAL<ri,t/T)dt  <3.12a> 

,^0 

No:*!  coefficient  a9  becomes  the  mean  of  the  function, 
since  WAL<0,t/'T)  »  1  over  the  normalized  interval  0it/TSl. 

The  Walsh  Series  can  be  defined  over  the  same  interval 
<0,T)  using  the  CAL  and  SAL  notations,  and  an  infinite  number 
of  terms:  t26:  p.2321 


<3. 15) 


1 


A  couple  of  simple  examples  follows  to  illustrate  the 
Mai sh  Series  expansion. 

a.  Expansion  of  8in<t>  in  a  Hal sh  Series 

The  sine  function  will  be  defined  over  one  periods 
with  a  period  of  2f.f  and  the  Hal  sh  functions  will  be  defined 
as  t/T  goes  from  0  to  1 . 


<3. 12a)  s 


Calculating  the  first  coefficient,  using  Equation 


aj  *•  1/21C 


a0  «■  l/2it 


*0*0 


<t)MAL<0,t/T)dt 


^2lC 

^sin<t)MAL 

p2lt 

^sin<t)  <  1) 


This  is  a  confident  calculation,  as  most  students 
Know  that  the  average  of  a  sinusoid  is  0! 

The  second  coefficient, 


ai  =  t/2X 


r2* 

jr 


<t)l*^t<t,t/T)dt 


aj  »  1/21C  •  C^sin<t><-’l)dt  +  \sin<t)<l)dt  ] 
aj  -  -2/t  «  -.43411 

A  lot  of  menial  integration  will  yield  only  these 

coefficients  of  any  practical  valuei 

ag  -  .24348 
ajg  ■■  .  12493 

Coefficients  numbered  2,4,9,  and  18  are  present 
but  have  values  less  than  0.04. 


Th«  -function  sin<t)  is  then  approximated  by 
sin<t>  a  ail<«^U<l,t)  +  a^L<5,t>  +  ai3MAL(13,t) 
where  the  coe-fficient  values  are  given  above. 

Figure  3.4  C23:  p.  ISl  shows  the  three  Mai sh 
functions  used  and  their  sum.  The  negative  sign  of  the 
coefficient  aj^^  has  been  used  to  invert  the  Ma1sh<l,t) 
function. 


Lwac(3,i) 


mL(l3,l) 


LT 


Figure  3.4. 


Sln(t>  with  a  Limited  Malsh  Series 


O-f  COUPS*,  th*  addition  o-f  more  o-f  the  terms  o-f 
the  Series  would  yield  a  better  approximation.  One  can  see 
that  representing  a  smooth  curve  using  rectangular  shapes 
would  require  a  large  number  o-f  terms  o-f  the  Walsh  Series. 

The  Walsh  Series  will  give  better  results 
representing  a  rectangular  function. 

b.  Expansion  of  a  Rectangular  Function 

A  rectangular  pulse  function  should  lend  itself 
nicely  to  representation  by  other  rectangular  functions. 

Indeed,  it  is  hoped  that  the  representation  might 
be  exact  or  with  small  error  in  only  a  few  terms. 

Consider  this  function: 

f<t)  a  1  0it<.2S 

8  .25<t<.73 

1  .75<tSi.8 


and. 


Proceeding  with  the  coefficients, 
f  <t)l«^L<0,t>dt 


--fi' 


25 

<  1>  <  l>dt 


1.0 

<  1)  <  Ddt 
.75 


-  0.5 


*2 


25 

<.  1)  <-l)dt 

8 


1.0 

( 1)<-I>dt 
75 


.5 


Coefficients  ap  a3,  and  above  will  all 
calculate  to  be  zero. 
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This  rectangular  time  function  *  which  could 
represent  two  radar  pulses,  can  be  represented  exactly  by  two 
Malsh  functions. 

f<t)  «  0.5*ML<O,t>  -  e.SMAL<2,t> 

Figure  3.S  shows  the  function,  the  Malsh 
functions,  and  the  sum  of  the  Malsh  functions.  This  is  quite  a 
difference  from  the  representation  by  Fourier  Series,  which 
would  require  a  lot  of  terms  to  present  an  inexact 
representation. 

c.  Mave-form  Synthesis 

These  results  obtained  above  indicate  that  a 
continuous  waveform  is  more  suited  for  Fourier  transformation 
and  a  discontinuous  waveform,  certainly  a  rectangular  one,  is 
more  suited  for  Malsh  transformat^ion .  This  conclusion  should 
be  easy  to  accept  with  examination  of  the  structure  of  the 
orthogonal  functions  <sine~cosine  and  Malsh). 

Figure  3.6  C23i  p.  333  shows  a  reconstruction  of 
two  waveforms  and  the  number  of  Fourier  and  Malsh  terms  csed. 
The  reconstruction  supports  the  conclusion  very  well.  In 
addition,  Beauchamp  £23i  p.39]  reports  that  the  same 
conclusion  can  be  drawn  when  the  Fourier  and  Malsh  Transforir; 
are  used  to  reconstruct  a  continuous  and  a  step  waveform. 

In  coming  to  this  conclusion,  Beauchamp 
considered  the  eff^c*  of  the  number  of  allowed  levels  for 
quantizing  the  waveform  and  the  sampling  interval.  C23:  p.351 
The  sampling  interval  determines  the  number  of  coefficients 
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Figure  3.4.  WAV«-form  Synthesis  Using  Walsh  Functions 


produced,  and  the  highest  Walsh  coe-fficient  -found  from  the 
Transform  determines  the  number  of  data  points  available  for 
any  quantized  level.  This  is  a  built-in  limitation  on  the 
accuracy  of  the  Walsh  Transform  that  the  Fourier  Transform 
doesn't  have. 

d.  Digital  Sampling 

A  theorem  similar  to  the  Sampling  Theorem  for 
sine-cosine  functions  is  applicable  to  the  Walsh  functions. 

The  minimum  sampling  rate  is  C23t  p.37] 

f^  =»  2^*^  (3. 16) 

with  k  being  the  power  of  2  that  represents  the  sequency 
bandwidth.  The  sequency  bandwidth  would  be  the  sequency  of 
the  highest  component  desired  to  be  represented. 

Suppose  a  CALC  12, t>  component  had  been  determined 
and  it  was  desired  to  include  this  component  in  the  bandwidth. 
The  sequency  of  this  component  is  12,  and  the  closest  power  of 
2  that  is  at  least  12  is  4  (2^  »  16).  Thus  K  -  4  and 

f,  *  *  32 

For  comparison,  a  cosine  of  frequency  12  would 
require  a  sampling  rate  of  24  <2Xfrequency> .  Aliasing  is 
equally  applicable  to  the  Walsh  functions  and  should  be 
considered.  (231  p.37] 

4.  The  Walsh  Transform 

Section  3  sfated  that  a  continuous  function  over  the 
interval  0<tSl  could  be  represented  by  a  sum  of  Walsh 
functions,  shown  in  Equation  3.11. 
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n»8 


<3. 11) 


wh«r« 


-Ji 


<t)MAL(n,t)dt 


<3.12) 


ThuSf  a  tr«ns-fortn  p«ir  can  b»  d»<fin*dt  C23:  p.403 


f<t)  «  Z  F<K)l«^L<K,t) 
n*0 


F<k)  -f'f<t)|oWL<k,t)dt 
J9 


<3.  17) 


<3.18) 


F<k)  could  b*  written •  -for  an  intorvat  BltST, 


F<k)  »  1/T  \  •f<t)l«WL<k,t)dt 


<3. 19) 


Now  «uppos«  th«  intorval  T  i«  divided  into 
N  parts,  letting  t  «  i4,  where  ^  is  a  small  time  increment  and 
i  is  an  integer.  The  aim  is  to  convert  the  integral  into  a 
summation  -for  discrete  computation. 


Then , 


F<k>  -  1/T 


l^L<k,t>dt 


with 


«  1/Nil  >  Z  •f<iiDMAL<k,iiD 


d  >■  dt 


N4  »  T 


similar  properties  as  the  discrete  Fourier  Coe-f-f icients  -from 


ev«n  or  odd  -functions.  That  is,  an  even  sequence  < symmetr ical  \ 

i 

around  its  midpoint)  will  transform  into  only  CAL  function 
coefficients,  and  an  odd  sequence  will  transform  into  only  SAL 
coefficients.  C23:  p.41] 

The  Discrete  Malsh  Transform  (DWT)  has  some  advantages 
over  the  Discrete  Fourier  Transform  <DFT)  because  it  involves 
only  additions  and  subtractions,  not  multiplications  as  in  the 
DFT <  Also,  the  transform  is  not  noisy  because  precise 
representations  in  the  digital  computer  .  ibie.  Precise 

representations  of  sine-cosine  functions  i  ■  OFT  is  not. 

Finally,  the  Walsh  transform  kernel  is  ±1  anv  is  considerably  ' 
easier  to  calculate.  £23:  p.42] 

An  important  point  of  the  transform  is  that  the  DWT  is  | 

its  own  inverse,  and  a  separate  inverse  transform  definition  ] 

j 

is  not  required.  See  Equations  <3.23>  and  <3. 24) .  I 

i 

a.  The  Fast  Walsh  Transform  j 

Computation  of  the  DWT  would  involve  j 

i 

mathematical  operations  of  either  addition  or  subtraction. 

i 

This  number  of  operations  has  been  decreased  by  using  the 
redundancy  in  the  Walsh  matrix  representation  of  the  Transform 
to  N1og2N  operations,  the  same  number  as  in  the  DFT.  £23: 
p.54]  Remember  that  the  DWT  involves  additions/subtractions, 
not  complex  multiplications/additions  as  in  the  DFT. 

The  Fast  Walsh  Transform  can  be  described  in  a 
signal  flow  diagram  that  results  in  a  "butterfly"  appearance 
just  as  in  the  FFT .  See  Figure  3.7  on  the  next  page.  Unique 
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*  •  * 
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m 


■fsAtures  of  th®  butterfly  in  the  final  stages  of  calculations 
offer  an  opportunity  to  calculate  the  FWT  *in-place“,  without 
using  computer  working  space.  t23:  p.SS]  The  result  is  not  in 
sequency  order,  however,  and  a  sorting  routine  must  be 
included  if  sequency  order  is  desired. 

/Appendix  A  gives  listings  of  the  FWT  programs  used 
in  this  thesis. 

b.  Malsh  Coefficients  and  the  Malsh  Matrix 

The  computation  of  the  coefficients  from  the  FWT 
is  very  fast  and  takes  advantage  of  the  properties  of  a 
Hadamard  matrix  [23:  p.S23 

The  Walsh  matrix  is  related  to  the  Hadamard  matrix 
and  can  be  used  to  determine  the  coefficients  of  the  Walsh 
Series.  C27:  p.41  Consider  this  matrix  equation, 

tanJ  =  1/N  »  CW]  •  Cxil  <3.27) 

where  IW]  is  a  N  X  N  matrix  of  the  Walsh  functions,  Lx^l  is 
a  column  vector  of  N  sampled  values  of  the  function,  and 
Ca,^],  n  >  0  to  N-i,  is  a  column  vector  of  the  Walsh 
coefficients. 

The  Walsh  matrix  can  be  thought  of  as  the  matrix 
of  sampled  Walsh  functions,  and  thus,  is  composed  of  plus  and 
minus  ones.  To  compute  the  coefficients,  simply  perform 
addition  and  subtraction  of  the  sampled  values  according  to 
the  row  of  the  Walsh  matrix  being  used.  An  example  will  show 
this  computation. 
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Let  the  number  o-f  sampler  in  the  column  vector  =  8 
*  N.  Writing  the  matrix  equation,  omitting  the  1/N  term, 


Can] 


CWl 


Cxj] 


<3.  i6> 


X0 

■ 

*1 

^2 

• 

*3 

X4 

• 

’‘S 

ll 

^7 

»0 

»2 

«4 

ag 

For  comparison,  use  N  sampled  values  of  the 
rectangular  function  on  paoe  73.  Thus, 

C  1  1  0  e  8  e  1  1  ] 

1  ^  0,1, 2, 3, Ha*, 7 

Now  evaluate  the  constants  af^  using  mat.  ix 
multiplication.  A  plus  sign  in  the  Walsh  matrix  equals  ^1, 
and  a  minus  si^ 
constant  by  N. 


Cxi]T 


means  - 1 . 

Don't  forget  to 

a0  ^  0.9 

a4  >  0 

ai  -  « 

*5  “  0 

a2  •  .9 

*6  ■  ® 

*3  -  • 

a7  -  0 

Comparing  with  the  value  of  9^2  9iven  in  the 
Walsh  Series  expansion  of  this  same  function  (page  73> ,  we 
note  a  sign  change.  The  reason  lies  in  the  use  of  positive 
phasing  for  the  Walsh  matrix  and  the  use  of  Harmuth  phasing  in 
the  Series  expansion  computation. 

The  same  answer  has  been  reached,  although  this  is 
the  discrete  case. 
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Matrices  can  also  be  used  to  show  the  Inverse 


Discrete  Malsh  Trans-form.  From  Equation  <3. 24)  , 

N-Cxj].  =  CWI-tanI  <3.28) 

»  CWI'CMI'Cx^I 

and 

Cxjl  »  1/N-CW3  .CW3 -CXil  <3.29) 

Equation  <3.29)  says  to  get  the  sampled  values  o-f 
the  original  within  a  constant  just  multiply  the 
coe-f-f icients  by  the  Malsh  matrix  CM]. 

c.  Input  Time  Shi-fts  and  Their  Ef-fect 

The  OFT  magnitude  is  invariant  to  the  phase  o-f  the 
input  signal.  Although  Malsh  Trans-form  signals  con-form  to 
Parseval-'s  Theorem,  the  DMT  is  not  invariant  and  the  same 
spectral  representation  cannot  be  achieved  independently  o-f 
the  phase  or  time  shi-ft  o-f  the  signal  .  C23:  p.42] 

The  DMT  ii  clme  invariant  when  the  time  shift  is 
obtained  by  dyadic  translation  C23:  p.43].  Namely,  i-f  x^, 
the  input  number  series,  is  shi-fted  to  where 

i^  a  madu1o-2  addition 


then , 

Zc*<k)  +  Z*a<k)  a  Xc*<k)  Xs*<k) 
i-f  the  coef-f icients  are  expressed  in  CAL  notation  <2^$  X^.) 
and  SAL  notation  <Z^,  X^)  .  See  Appendix  C  -for  modu1o-2 
addl tion . 


The  importance  o-f  variations  with  time  shi-ft  can 
be  lessened  by  considering  sums  o-f  squares  o-f  coe-f-f  icients. 
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This  is  one  definition  of  the  Walsh  Power  Spectrum,  and  is 
considered  in  the  next  section. 

5.  The  Walsh  Power  Spectrum 

Spectral  analysis  using  sequency  has  some 
advantages  over  analysis  in  the  -frequency  (Fourier)  domain. 

A  defined  power  spectrum  from  the  Walsh  functions  is  easily 
applied  to  discontinuous  or  time  limited  functions,  for 
example,  and  it  is  possible  to  have  a  sequency  limited 
spectrum  for  a  time  limited  function  in  contrast  to  the 
Fourier  Spectrum.  t23t  p.871 

Although  there  are  several  derivations  of  the 
Walsh  power  spectrum  C23:  p.91],  a  more  normal,  i.e. 
analogous  to  Fourier,  approach  will  be  used  to  make  a  spectrum 
definition.  As  there  are  readily  available  methods  (FVfT)  for 
calculating  the  Walsh  coefficients,  a  combination, 
specifically  the  sum  of  the  squares  of  related  coefficients, 
would  give  an  easy  to  use  definition.  Also,  as  stated 
earlier,  by  using  the  sum  of  the  squares  of  certain 
coefficients,  the  time  shift  variance  that  plagues  the  Walsh 
Transform  can  be  lessened. 

Define,  then,  the  Power  Spectral  Density  <PSD> 
coefficients  of  the  Walsh  spect'^um.  123:  p.l00] 

P<0)  -  Xc*<0,t>  <3.31) 

P<k>  «  Xc*<k)  +  X^a<k) 

P<N/2)  »  X^a<N/2) 

k  «  1,2, . . . ,<N/2  -  1) 


Not*  that  <N/2  +  1)  spectral  points  are 
generated.  This  method  is  analogous  to  Fourier  power  spectral 
analysis,  where  the  power  coefficient  is  the  f>um  of  the 
squares  of  the  real  and  imaginary  parts  of  the  complex  Fourier 
Transform  coefficients. 

Beauchamp  C23:  p.i80]  states  that  taking  the 
square  root  of  P(k>  does  not  give  a  sequency  amplitude 
spectrum,  but  other  references  C26:  p.233]  use  as  the 
definition  of  the  sequency  spectrum. 

Another  spectrum  can  be  defined  using  the  Malsh 
Coefficients.  The  name  Group  Spectrum  is  the  author'' s  term 
and  will  be  used  here.  Beauchamp  £23:  p. 1063  refers  to  it  as 
"the  odd-harmonic  sequency  spectrum",  and  Campanel 1 a/Robinson 
t26t  p.234]  refer  to  this  spectrum  definition  as  a  generalized 
Wal sh-Four ier  spectrum  C2B]. 

The  power  content  of  individual  groups  of 
frequency  components  has  been  shown  to  be  equal  to  groups  of 
components  in  the  sequency  domain.  Any  group  of  components, 
in  either  domain,  is  a  collection  of  the  orthogonal  components 
that  make  up  the  original  signal  or  function. 

Call  groups  of  sequency  components  6n f  with  £nl 
an  integer,  and  consider  a  discrete  signal  made  up  of  a 
sequence  Csk>  of  length  N. 

The  total  energy  of  the  signal  is,  by  Parseval's 
Theorem,  equal  to  the  sum  of  the  orthogonal  components  t2‘,  ]. 
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<Th*sc»  orthogonail  s»qu«ncy  componsnts  havtt  bean  grouped  into 
Gf,.)  Then, 

N-l  m 

l/N  •  »  rOp  <3.32) 

k“0 '  n“0 

and  the  number  o-f  groups  Cml  is  given  by 

m  *■  I  +■  1og2N  <3.33) 

The  groups  69  and  6<1og2^)  represent  the  power 
content  o-f  the  d.c.  component  and  the  -folding  -frequency 
component,  respectively.  The  other  Groups  represent  the  power 
content  o-f  the  component  atn^  1,  n^2,  n«3,  etc.,  plus 
all  odd  harmonics. 

G9  »  P*<0)  <3.34) 

Gi  -  pa<l)  +  P*<3>  ♦  P*<5)  +  ... 

G2  -  P2<2)  -»•  pa<d>  *■  Pa<10>  +  ... 

Qg  »  pa<4)  +  pa<i2)  -••  pat<20)  ■••... 

G4  »  Pa<8)  ♦  pa<24>  *  P2<40)  +  ... 
and  so  on,  until 

G<1og^>  •  P2<N/2> 

The  plot  of  G,)  versus  n  gives  a  discrete  spectum 
that  is  unique  for  a  certain  time  sequence  C26:  p.234],  and 
invariant  to  a  time  shift  of  the  input  C23:  p.l06). 

This  spectrum  is  highly  compressed  <m  points 
instead  of  N/2  -*■  1  points),  and  therefore  doesn^t  present  all 
characteristics  of  the  data  as  well  as  the  previously  defined 
spectrum.  C23:  p.l061 


The  application  o-f  the  Walsh  Transform  to 
continuous  or  discontinuous  functions  is  easily  accomplished 
using  the  FWT  programs  provided  in  Appendix  A.  The  programs 
input  a  number  series  <which  could  represent  a  sampled 
waveform)  and  compute  the  FWT.  With  these  coefficients,  the 
Power  Spectral  Density  and  the  Group  spectrum  can  then  be 
calculated.  The  FWT  will  be  applied  to  the  previous  examples 
and  the  results  will  be  compared  and  analyzed  in  terms  of  the 
elaborated  theory  in  this  chapter. 

a.  The  R4T  of  a  Sinusoid. 

Figure  3.8  gives  a  normalized  graph  of  the 
coefficients  of  the  FWT  of  one  cycle  of  a  sinusoid  in  the 
Walsh  interval  0  to  i .  The  actual  magnitudes  can  be 
determined  from  the  scale  listed  on  the  vertical  axis. 

Running  either  program  in  Appendix  A  will  give  a  listing  of 
the  coefficient  magnitudes. 

Note  that  the  maximum  coefficient  (normalized 
to  1)  is  a  WAL(n,t)  •  MAL<l,t>.  Equation  (3.4)  of  the 
chapter  reveals  there  is  a  SAL<lyt>  with  a  sequency  of  1  at 
this  Walsh  index.  Thus,  the  one  cycle  sin  wave  in  the  Walsh 
interval  0  to  1  has  a  distinguishing  feature  at  a  sequency  of 
1.  Perhaps  not  surprising  since  i*  is  Known  that  frequency  is 
»  special  :r:%.aisure  of  sequency,  and  the  one  cycle  in  this 
interv^.l  can  thought  of  as  having  a  frequency  of  1. 


V  Al«o  not*  the  other  prominent  components  at 
n  »  9,  9,  13f  and  29,  with  sequencies  of  3,  5,  7,  and  IS, 
respecti,*vet y.  Our  previous  calculations  with  the  Malsh  series 
(page  71)  predicted  conspicuous  coefficients  at  n  »  5  and 

:  13.  V  '  '■ 

V  (3.31)  to  these 

coefficients  results  in  Figure  3.9,  the  Power  Spectral  Density 
(PSD).  The  )Aax i mum  value  of  power  is  at  a  sequency  of  1,  and 
there  is  no  DC  power  term  (sequency  of  0> .  Harmonic  power  is 
located  at  sequencies  of  3,  7,  and  15. 

These  odd  number  power  components  combine 
into  f  Single  component,  Gj,  in  the  Group  spectrum.  Figure 

*  t  I  4 

3.10.  See  Equation  <3.34> . 

An  increase  of  frequency  gives  similar 
results.  Figures  3.11  to  3.13  show  the  FWT,  PSD,  and  Group 
Spectrum , of  three  cycles  of  a  sinusoid  in  the  Malsh  interval  8 
to  1.  The  maximum  FMT  coefficient  lies  at  n  »  5  Ca  SAL(3,t> 
component  with  a  sequency  of  33,  with  other  major  components 
at  n  »  1,  9,  13,  21,  25,  29,  53,  and  57.  All  of  these 
components  are  at  odd  sequencies,  so  the  power  will  be  located 
in  the  odd  number  sequencies  of  the  PSD,  and  the  odd  numbered 
components  will  combine  into  one  component  of  the  Group 
Spectrum.  Again,  see  Equations  (3.34) . 

The  FMT,  PSD,  and  Group  Spectrum  of  eight 
cycles  of  a  sine  wave  in  the  Malsh  interval  is  analyzed  in 
Figures  3.14  to  3.16.  The  FWT  is  considerably  simpler,  with 
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Figur*  3.10.  Group  Spectrum  o-f  1  Cycle  o-f  4i  Sine  Function, 
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Figure  3.14.  FWT  o-f  8  Cycles  of  a  Sine  Function. 
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only  4  components  at  PSD  sequencies  o-f  8  and  24.  The  PSD 
shoM4s  the  peak  power  at  a  sequency  of  8  and  a  notable  harmonic 
at  24,  The  Group  Spectrum  -follows  from  Equations  3.34. 

A  few  conclusions  about  the  transform  of 
a  sinusoid  could  be  drawn.  The  peak  coefficient  of  the  FWT 
falls  at  a  Walsh  index  that  translates  to  a  sequency  that  is 
equal  to  the  sinusoidal  frequency.  The  PSD  then  has  a 
component  that  is  maximum  at  this  sequency.  Irt  addition,  this 
maximum  power  component  is  usually  accompanied  by  notable 
components  at  harmonics  of  this  sequency. 

The  Group  Spectrum,  although  highly 
compressed  and  uncomplicated,  requires  a  knowledge  of  the  PSD 
for  it  to  be  useful.  Note  that  the  Group  Spectrum  of  the  1 
cycle  sine  wave  is  equal  to  the  the  Group  Spectrum  of  the  3 
cycle  wave.  An  educated  guess  says  that  the  Group  Spectrum 
for  all  odd  numbers  of  cycles  would  be  the  same,  that  of  G^ , 
because  only  odd  numbered  components  CEq.  3.34)  of  sequency 
are  present.  Each  component  of  the  Group  Spectrum  could 
represent  a  series  of  sinusoidal  frequencies,  with  no  way  to 
distinguish  the  actual  frequency  of  the  sinusoid  without  prior 
knowledge  of  the  PSD. 

b.  The  FWT  of  a  Rectangular  Function. 

Consider  the  FWT,  PSD,  and  Group  Spectrum  of  a 
4  hertz  square  wave.  Knowing  that  the  Walsh  function  lends 
itself  readily  to  the  reproduction  of  such  a  wave,  simple 


9? 


results  -for  graphs  o-f  the  transform  and  related  computations 
would  be  expected. 

Indeed  this  is  the  case.  Figures  3.17  and 
3.18  show  the  simplicity  o-f  representation  -for  rectangular 
■functions.  The  FWT  component  at  n  =  7,  sequency  *  4,  is  the 
only  coe-f-f icient  in  Figure  3.17.  (Note  that  a  4  hertz  square 
wave  is  WAL(7,t)).  All  of  the  power  lies  at  a  sequency  of  4 
(Figure  3.13),  and  in  a  (sroup  Number  3  (Eqn.  3.34).  With  a 
proper  number  of  samples,  any  square  wave  of  a  particular 
frequency  would  be  represented  with  similar  characteristics. 

A  rectangular  function,  page  73,  composed  of 
square  pulses  was  expanded  in  a  Walsh  Series  in  this  chapter. 
Consider  the  FWT,  PSD,  and  Group  Spectrum  of  the  function. 

Figures  3.19,  3.20  and  3.21  show  the  results 
of  the  computation.  These  graphs  are  not  normalized  in 
order  to  show  the  similarities  with  the  previously  computed 
example  on  page  73.  The  FWT  shows  2  coefficients  of  0.5 
magnitude,  at  n  »  8  and  n  ™  2,  matching  the  results 
determined  before.  Since  the  FWT  outputs  coefficients  based 
on  positive  phasing,  this  result  is  correct. 

The  power  contained  in  the  function  is  in  the 
DC  term  and  the  first  sequency  component  as  shown  in  Figure 
3.20.  These  sequency  components  result  in  2  components  in  the 
Group  Spectrum.  See  Equations  3.34. 
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Figure  3.21.  Group  Spectrum  o-f  2  Rectangular  Pulses. 
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c.  The  FWT  and  a  Input  Time  Shi+t. 

Section  4-c  stated  that  the  Discrete  Walsh 
Transform  is  not  invariant  to  a  time  shift  or  phase  shift  of 
the  input.  Let  sin<t)  be  circularly  shifted  i£/4  to 
sin<t-lC/4).  Now  examine  Figures  3.22,  3.23.  and  3,24  and 
compare  with  Figures  3.8,  3.9,  and  3.10. 

Note  the  shifted  FVrr  coefficients  pattern  is 
more  complex.  A  sizable  component  now  lies  at  n  ~  2.  This  is 
a  CAL<l,t>  component  with  a  sequency  of  1.  The  unshifted  FWT 
graph  contained  only  a  negligible  CAL(l,t>  component. 

Double  components  at  n  5  and  6,  13  and  14,  and  29  and  38 
represent  sequencies  of  3,  7,  and  15,  r espec t i vel y ,  the  same 
components  that  wore  contained  in  the  unshifted  sinusoid,  but 
the  magnitudes  of  the  components  at  the  same  index  CnJ  are  not 
equal . 

Although  the  FWT  magnitudes  are  different  at 
each  index  Cn3  of  the  shifted  and  unshifted  sinusoids, 
remember  that  adjacent  components  combine  to  form  a  PSD 
coefficient  of  a  particular-  sequency.  The  components  vary  in 
a  reciprocal  manner,  so  when  they  are  squared  and  summed,  the 
time  variant  effect  is  indeed  lessened  1235  p.891.  The  result 
is  a  PSD  that  has  sequency  components  that  are  very  close  to 
being  time  invariant. 

One  expectation  that  was  revealed  to  be  true 
was  the  time  invariance  of  the  Group  Spectrum.  Beauchamp  C23; 
p.l96]  states  this  spectrum  is  time  invariant  to  time  shifts 
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Figure  3.22.  FWT  of  a  Shifted  Sinusoid. 
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Figure  3.23.  PSD  o-f  a  Shifted  Sinusoid. 
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Figure  3.24.  Group  Spectrum  o-f  a  Shifted  Sinusoid. 
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o-f  the  input,  and  examination  o-f  Figure  3.23  and  its  magnitude 
gives  an  af+irmative  response  to  his  assertion. 

The  rectangular  pulse  function  will  positively 
demonstrate  the  effects  of  time  shift.  The  FWT,  PSD,  and 
Group  Spectrum  of  Figure  <3.2Sa)  is  shown  in  Figures  3.19, 
3.20,  and  3.21.  The  FWT  and  associated  spec tr urns  of  Figure 
<3.25b>  is  shown  in  3„26,  3.27,  and  3.28. 

Quite  a  difference  in  the  FWT  of  the  shifted 
function  can  be  noticed.  The  coefficient  pattern  is  much  more 
‘spread  out",  and  the  magnitudes  of  the  prominent  coefficients 
at  n  as  0  and  1  are  less  in  the  shifted  rectangular  function. 

The  PSD's  arc  different  also.  With  the  power 
being  spread  among  more  sequencies  in  the  shifted  rectangular 
function,  the  magnitudes  of  the  power  coefficients  at  0  and  1 
are  reduced. 

The  Group  Spectrums  are  not  alike  in  this 
case.  The  two  input  sequences  to  the  FWT  program  are  made  up 
of  samples  of  the  original  function  and  the  shifted  function. 
These  sequences  are  not  al iKe  and  generate  different  Group 
Spectrums. 

Has  Beauchamp‘'s  assertion  [23:  p.l061  been 
violated?  The  answer  is  no.  In  the  sinusoidal  shift,  it  was 
assumed  that  the  function  was  periodic  outside  the  interval 
0<t<l  over  which  the  Walsh  functions  are  defined  and  over 
which  the  function  was  taken.  When  the  sinusoid  was  shifted, 
the  part  shifted  out  of  the  interval  "wrapped  around"  into  the 


Figgr®  3.25.  Two  Rectangular  Pulses,  Shifted  and  Unshifted 
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Figur#  3.27.  PSD  of  Shifted  Rectangular  Pulses. 
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■first  part  o-f  the  interval  (circular  shi-ft)  ,  resulting  in  the 
shi-fted  -function  having  the  same  power  as  the  unshifted 
function . 

However,  the  rectangular  function  wasn^'t 
defined  as  a  periodic  function,  merely  two  pulses  in  an 
interval.  When  the  rectangular  function  was  shifted,  the  part 
that  shifted  outside  the  interval  wasn't  wrapped  around  to  the 
front  (not  circularly  shifted)  and  was  lost,  decreasing  the 
total  power  found  in  the  shifted  function.  The  PSD  or  the 
Group  Spectrum  could  not  be  the  same  as  the  unshifted 
function,  since  the  amount  of  total  power  is  different  in  each 
of  the  functions. 

A  conclusion  is  made  that  the  Group  Spectrum 
is  time  invariant  to  phase  shifts  in  the  input  for  periodic 
circular  shifted  functions  only. 


IV.  THE  UWLSH  TWWSFORM  DEI^fTERLEAVIN6 


Th«  -first  two  chapters  presented  some  background  -for  the 
appreciation  o-f  the  role  of  ESM  and  the  basic  components  o-f  a 
typical  ESM  system.  Understanding  the  tools  and  purpose  of  the 
trade,  though,  only  gives  a  small  insight  into  the  nature  of 
probi em. 

There^s  a  lot  of  electromagnetic  radiation  out  there!  The 
combat  environment  presents  an  overwhelming  amount  of  energy 
that  has  to  be  collected,  analyzed,  and  sorted.  The  data 
stream  pours  into  the  system  from  all  angles,  with  numerous 
frequencies  and  amplitudes.  Corruptions,  disturbances, 
reflections,  and  missing  pulses  make  the  deinterleaving 
process  a  most  difficult  problem.  Today's  techniques, 
described  in  Chapter  Two,  do  a  good  Job,  but  improvements  are 
needed  and  sought  after. 

Unique  in  many  respects,  the  Walsh  Functions,  which  are 
described  in  Chapter  Three,  are  examined  for  their  usefulness 
in  deinterleaving  pulse  trains  into  separate  chains  of  pulses 
with  different  pulse  repetition  intervals.  More  specifically, 
this  chapter  attempts  to  determine  if  there  is  a  PRI/PRF 
recognition  feature  in  the  Malsh  Transform  of  the  dat?. 
presented  by  an  ESM  receiver. 


A.  REPRESENTATIC3N  OF  ESM  RECEIVER  DATA 


The  operating  element  o-f  the  typical  deinterleaving 
process  is  a  computer  data  word  <or  words)  o-f  some  format 
appropriate  to  the  system. 

The  data  word  contains  values  of  the  initial  parameters 
measured  by  the  system,  and  basically  these  values  are 
examined  for  matches  in  previously  analyzed  pulses.  Software 
processing  systems  make  these  comparisons  and  places  words  of 
similar  parameters  into  cells  or  histograms,  patiently  waiting 
until  enough  pulse  information  can  be  found  to  make  an 
identification.  Other  algorithms  process  the  data  word  and 
make  calculations  of  PRI's,  adding  a  new  sorting  parameter 
that  is  used  in  the  deinterleaving  process. 

The  investigation  of  deinterleaving  with  the  FWT  will  need 
a  different  symbology  for  the  receiver  data,  since  the  FWT 
operates  on  a  series  of  N  numbers,  not  a  data  word,  to  obtain 
the  transform. 

Imagine  an  absolute  time  line  that  consists  of  some  point 
in  the  past,  and  ends  at  the  present.  Consider  a  received 
pulse  as  being  represented  by  a  vertical  line  stationed  at 
the  time  of  arrival  <TOA)  of  that  pulse.  The  height  of  the 
line  could  be  made  proportional  to  the  amplitude  of  the 
received  pulse.  The  point  on  the  time  line  where  the  pulse 
"stands*  is  the  time  of  arrival  of  the  pulse. 

Figure  4.1  looks  like  a  messy  situation,  but  there  are 
really  only  four  pulse  trains  of  different  PRI's  on  the  line. 
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Figure  4.1.  TO^  representation  o-f  received  pulses. 

PRI's  could  be  calculated  by  subtracting  successive  TOA's 
o-f  pulses  in  the  same  pulse  train. 

PRI  1  =  T0Aj2  -  TOAii 
PRI2  =•  TOA22  ”  TOA21 
PRI 3  =  TOA32  -  T0A3i 
PRI 4  =»  TOA42  *■  TOA41 

This  procedure  is  used  uihen  other  parameters  are  present 
to  identify  the  pulses. 

Keep  in  mind  that  there  is  a  steady  stream  of  these  pulses 
into  the  ESM  system.  In  order  to  maKe  things  manageable,  the 
number  of  pulses  being  examined  has  to  be  reduced. 

Lemley  129  i  p.l2]  calls  this  reduction  mechanism  the 
■presorting  aperture.*  It  is  of  variable  width  to  accomodate 
high  and  low  density  environments  and  is  similar  to  the  time 
slice  concept  discussed  under  deinterleaving  in  Chapter  Two. 

It  reduces  the  number  of  pulses  being  processed  at  one  time. 

Another  way  of  thinking  of  this  reduction  concept  is  to 
visualize  a  “windc«»*“  on  the  stream  of  pulse  lines.  Only  a 
certain  number  of  pulse  lines  can  be  present  in  the  window  at 
anytime,  and  deinterleaving  operations  are  carried  out  on 


118 


these  pulse  Ines.  The  window  can  be  divided  into  small  time 
intervals,  which  could  be  analogous  to  the  resolution 

specification  o-f  a  receiver  (usually  100-200  nanoseconds).  If 
a  pulse  line  is  present  in  the  time  interval,  a  number 
proportional  to  the  amplitude  is  generated  for  the  pulse  line 
representation.  If  not,  a  zero  is  used. 

The  author  prefers  to  visualize  that  the  receipt  of  a 
pulse  triggers  the  collection  of  the  next  N  time  intervals  for 
analysis  and  deinterl eaving.  Each  interval  could  contain  a 
pulse  or  not. 

B.  CONSTANT  fiWPLlTUDE  TOA  STRINGS  THE  FWT 

This  chapter  will  work  mostly  with  N  =*  64  intervals  and 
pulse  line  representations.  N  must  be  a  power  of  two  to  be 
used  in  the  FWT, 

The  pulse  line  representation  can  be  symbolized  for 
computer  work  with  the  BASIC  language  DATA  statement.  Using 
the  statement  on  the  64  pulse  lines  in  the  window  gives  a 
representation  that  accounts  for  the  two  variables  that  are 
present  in  the  pulse  line  stream,  the  amplitude  and  TOA  of  the 
pul ses. 

Such  a  DATA  line  for  Figure  4.1  would  be 

DAT  A  1 ,  .  3 ,  75 ,  .3,0,  .5,  .  75 , 1,0,  .3,  .5,  .3,0,  .3,  1,  .5,0,  .3,  .75 

<4.  1> 

Each  number  represents  the  amplitude  of  a  recognized  pulse 
that  occupies  a  TOA  time  interval . 
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Think  o-f  it  in  one  of  two  ways: 

1.  The  receiver  outputs  a  measured  value  of  a  pulse 
amplitude  each  TOA  resolution  time  interval.  If  no  pulse  was 
present,  it  outputs  a  zero.  Or, 

2.  The  receiver  outputs  a  pulse  amplitude  only  upon 
detection  and  analysis  of  a  pulse.  Time  intervals  of  TOA 
resolution  width  in  between  pulse  detections  can  be 
represented  with  zeroes. 

Either  way,  the  DATA  line  looks  like  (4.1). 

Assume  for  simulation's  sake  that  the  period  between  DATA 
numbers  is  a  known  value  of  time,  and  the  DATA  line  represents 
the  output  of  an  ESM  receiver  to  the  preprocessor.  The 
numbers  indicate  amplitude  of  the  received  pulses,  and  their 
distance  in  interv:^ls  between  similar  pulses  is  a  measure  of 
their  PRI's. 

Things  can  be  simplified  a  bit  further  by  considering  all 
pulse  lines  of  the  TOA  string  to  be  of  the  same  amplitude. 

Only  one  variable  would  then  be  present  in  the  string,  the 
TOA. 

A  DATA  line  representation  of  a  string  of  this  type  is 
DATA  1,1, 1,1, 9, 1,1, 1,8, I, 1,1, 0,1, 1,1, 0,1,1,  (4.2> 

In  short,  this  TOA  string  and  its  DATA  representation  is 
merely  an  indication  of  the  rec*'ipt  of  a  'tulse  or  not  during 
the  resolution  time  interval . 

A  number  of  Fast  Nal sh  Transforms  were  run  on  simulated 
DATA  TOA  strings,  and  plots  of  the  coefficients  were  made.  The 
Power  Spectral  Density  and  Group  Spectrum  coefficients  were 
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a1  calculated  and  plotted.  All  o-f  the  Group  Spectrum 
plots  are  not  included,  as  the  emphasis  o-f  this  examination 
was  placed  on  the  FMT  and  PSD.  Some  o-f  them  are  discussed  in 
Section  C  o-f  this  chapter. 

The  plots  are  grouped  by  type  <FWT  or  PSD)  in  Appendix 
D.  Each  plot  is  normalized  by  dividing  all  coe-f-f icients  by 
the  maximum  coefficient.  Actual  values  of  the  coefficients 
can  be  determined  with  the  scale  on  the  vertical  axis.  The 
maximum  component's  magnitude  is  shown  in  parenthesis. 

These  plots  show  the  index  number  and  sequency  of  the  FWT 
and  PSD  components,  respectively,  but  only  their  relative 
magnitude  with  respect  to  the  maximum  component.  Keep  that  in 
mind  when  examining  the  plots.  Where  important,  the  FVTT  and 
PSD  are  replotted  with  the  actual  magnitudes  on  the  vertical 
axis. 

These  DATA  representations  simulate  single  PRI  pulse 
trains  received  by  an  ESM  system.  Note  the  number  of 
intervals  between  successive  pulses  is  indicated  in  the  title 
of  the  plot.  To  calculate  a  simulated  PRI  of  the  train, 
simply  multiply  the  number  of  intervals  between  sucessive 
pulses  by  the  time  interval  resolution,  a  Known  value  and  a 
function  of  the  receiver. 

The  plots  are  examined  by  type  in  the  following  sections, 
and  conclusions  are  reached  and  reported.  These  single  PRI 
representations  are  the  building  blocks  for  the  interleaved 
pulses  that  are  investigated  later  in  the  chapter. 


1.  Fast  Mai  sh  Trans-forms  o-f  the  TOA  Strings 


An  examination  o-f  the  FWT-'s  separated  them  into  two 
groups  o-f  similar  characteristics. 

'*‘he  first  group  is  composed  of  FWT's  with  intervals 
between  pulses  that  are  even  numbers.  Examine  each  plot 
in  Appendix  D  (page  176  to  i97>  and  note  that  each  even 
interval  FWT  is  symmetr i cal  1 y  even  about  a  point  midway 
between  n  =  31  and  32.  Also,  if  the  even  number  is  a  power  of 
two  (p.o.t.),  the  FMT  is  considerably  less  complex.  Each 
successive  p.o.t.  interval  is  composed  of  the  addition  of  a 
pair  or  pair  of  coefficients  to  the  previous  p.o.t.  FV4T,  and 
its  magnitude  is  one  half  the  previous  p.o.t.  FWT. 

For  example,  one  pair  of  coefficients  is  added  to  TCA 
TAS  2  at  n  =*  31  and  32  (sequency  of  16)  to  form  TOA  TAG  4, 

Two  pairs  of  coefficients,  <n  =  15  and  16,  n  =  47  and  48)  are 
added  to  TOA  TAG  4  to  make  TOA  TAG  8,  4  pairs  to  make  TOA  TAG 
16,  etc.  At  N  =  64,  enough  pairs  of  coefficients  have  been 
added  to  have  one  coefficient  for  each  Cnl. 

The  odd  numbered  intervals  do  not  have  symmetrical 
characteristics.  In  addition,  the  highest  relative  magnitude 
Walsh  index  component  or  components  seem  to  bear  no  relation 
to  the  interval  number  and  thus  the  PRl  of  a  pulse  train.  (An 
unfortunate  conclusion  for  PRI  recognition.) 

Close  scrutiny  of  both  appearance  groupings  was  made 
in  an  effort  to  notice  a  feature  that  is  distinctive  to  the 
particular  simulated  PRI.  There  doesn't  seem  to  be  any 
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present,  either  in  m&gnitude  or  coe-f-f icient  distribution,  that 
is  obvious  without  a  detailed  examination. 

0-f  initial  interest  to  the  author  was  the  magnitude  o-f 
the  MAl.<9,t)  component.  This  magnitude  seemed  to  be  related 
to  the  number  o-f  resolution  intervals  between  the  pulses.  It 
holds  -for  most  but  not  all  cases.  Consider  Table  4.1  which 
lists  the  TAG  number  and  the  WAL(9,t>  component. 


Table  4.1.  Amplitude  o-f  MAL<8,t>  and  TAG  No.  Inverse 


1  TAG  no. 

WAL(9,t)  amplitude 

TAG  no.  inverse 

1 

1.8 

1.0 

2 

0.5 

0.5 

3 

0.34375 

8.3333 

4 

8.25 

0.25 

5 

8.283125 

0.20 

6 

0.171875 

0 . 1667 

7 

8 . 15625 

0 . 1429 

8 

8 . 125 

0 . 125 

9 

0.125 

0.1111 

18 

0 . 189375 

0 . 1 

11 

0.09375 

0 ,0909 

12 

0.09375 

8 .0833 

13 

0.978125 

0 .076923 

14 

0 .078125 

0.071429 

IS 

0 .078125 

0 .0667 

16 

0.0625 

0.0625 

19 

0.0625 

0.0526 

22 

0,046875 

0.045455 

25 

0.046875 

8  .84 

32 

0.83125 

0 .03125 

48 

0.03125 

8.025 

51 

0 .03125 

0.019608 

64 

0 .015625 

0 .0 15625 

For  instance,  TOA  TAG  2  has  two  intervals  between 
pulses,  and  the  magnitudeo-f  MAL(0,t>  is  8.5,  or  the  inverse  o-f 
2.  TOA  TAG  3  has  a  magnitude  o-f  .34375,  which  is  close  to  the 


inverse  of  3.  TOA  TAG  4  has  a  WAL<8,t>  amplitude  of  .25, 
which  is  one  divided  by  four,  the  number  of  intervals  between 
pulses  in  the  representation. 

If  there  was  a  relationship  between  the  number  of 
intervals  between  the  pulses  in  the  DATA  line  and  the 
amplitude,  then  perhaps  this  would  carry  over  into  the 
interleaved  pulse  trains  composed  of  these  single  PRI 
representations.  Exceptions  are  evident,  though,  such  as  TOA 
TAG  9,  which  has  the  same  magnitude  as  TOA  TAG  8,  although  ihe 
coefficient  distribution  is  not  the  same.  This  implies  that 
any  recognition  feature  will  probably  depend  upon  two 
variables,  perhaps  the  amplitude  and  the  sequency  of  the 
coefficients. 

2.  Power  Spectral  Densities  of  the  TOA  Strings 

An  examination  of  the  PSD's  in  Appendix  D,  page 
198  through  page  219,  will  yield  similar  conclusions. 

PSD's  of  even  intervaled  TOA  strings  are  symmetrically 
even  about  a  sequency  of  16.  PSD's  of  even  numbers  that  are 
also  a  p.o.t.  are  considerable  less  complex,  and  are  one 
fourth  the  magnitude  of  the  previous  p.o.t.  PSD  (except  for 
TOA  TAG  2  to  TOA  TAG  4) . 

The  odd  numbered  PSD's  have  no  symmetry,  and  sequency 
components  that  stand  out  in  relative  magnitude  do  not  bear  a 
relationship  to  the  simulated  PRI  of  the  string.  (Remember 
the  number  of  intervals  between  the  TAG'S  is  a  measure  of  the 
PRI  of  the  simulated  pulse  train,  when  the  time  interval 
between  numbers  is  Known.) 
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3>  General  Comments 


The  entire  purpose  o-f  this  examination  has  been  to 
recognize  similarities  or  relationships  between  succcessive 
TOA  TAG  representations  in  order  to  recognize  the  same 
-features  (i-f  they  exist)  in  an  interleaved  pulse  train  DATA 
line  representation. 

At  this  point,  the  results  are  very  general  and 
an  optimistic  mood  regarding  any  usable  -features  might  belong 
only  to  the  optimist.  It  seems  that  the  even  numbered  graphs 
<FWT  and  PSD  )  are  the  only  graphs  that  show  any  usable 
■features  that  could  be  used  -for  recognitor.  One  shouldn't 
give  up  all  hope,  though. 

The  author  believes  it  to  be  possible  to  develop  a 
computer  algorithm  that  would  recognize  and  determine  the  TOA 
TAG  number  (and  thus  a  PRI  o-f  a  simulated  pulse  train)  -for  a 
given  single  pulse  train  DATA  line  representation.  Enough 
-features  exist  in  the  FWT  and  PSD  that  could  be  used  to 
process  a  particular  single  PRI  representation  with  a  Known 
resolution  time  interval  and  output  the  correct  PRI.  The 
pr ogt  am  would  have  to  determine  the  PRI  by  matching 
characteristics  (sequency  components  and  magnitudes)  o-f  the 
FWT  and/or  PSD  to  a  data  bank  composed  o-f  characteristics  o-f 
the  FWT's  and  PSD-'s  of  single  pulse  train  representations. 

For  example,  it  could  first  take  into  account  the  symmetry  of 
the  FWT/PSD.  If  it  is  symmetrically  even,  then  it  isn't  an 
odd  interval  TOA  TAG. 
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The  algorithm  would  be  reasonably  complex  for  just 
this  simple  TOA  TAG  representation  of  a  single  pulse  train  of 
a  particular  PRI .  But  it  could  be  done. 

Since  it  addresses  only  single  PRI  TAG 
representations,  the  algorithm  would  only  be  effective  if  you 
could  convince  an  opponent  tu  use  only  one  radar  at  a  time 
against  the  ESM  platform.  Obviously,  this  isn''t  practical. 

More  realistically,  an  examination  of  interleaved 
pulse  trains  and  a  determination  of  a  PRI  recognition  feature, 
if  one  exists,  would  be  in  order.  This  will  be  done  in 
Section  D.  after  discussion  of  the  effects  of  a  time  shift  in 
the  TOA  TAG  string. 

Seme  limitations  of  the  DATA  line  symbology  should  be 
stated.  Using  N  »  ^4  and  a  resolution  between  numbers  of  200 
nanoseconds  means  the  window  length  is  less  than  2 
microseconds  long.  This  is  not  a  realistic  figure,  since  a 
high  density  of  10*  pulses  per  second  gives  an  average  of  1 
pulse  every  microsecond.  The  use  of  N  «  1024  for  an 
interleaved  DATA  line  representation  might  be  far  better  in 
the  numbers  of  coefficients  (and  more  realistic  in  length  of 
the  window),  but  unless  an  identifiable  PRI  feature  is 
present,  the  extra  data  is  exactly  that,  extra  data. 

An  additional  limitation  of  this  representation 
concerns  coincident  pulses.  Pulses  arriving  at  the  receiver 
at  the  same  time  (falling  in  the  same  resolution  time 
interval)  are  not  individually  represented  in  the  DATA  line 
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symbology.  For  this  reason ^  one  cannot  represent  a  TOA  TAG  2 
and  any  other  even  numbered  TAG.  This  limitation  is  a 
problem  that  will  a-f-fect  the  FWT  and  PSD  o-f  the  interleaved 
pulse  train  representation. 

C.  TIME  SHIFTED  TOA  STRINGS 

Chapter  Three  has  stated  that  the  FWT  and  PSD  are  not  time 
invariant  to  circular  shifts  of  the  input  number  series. 
Figures  4.2  through  4.3  are  FWT's  of  shifted  TOA  TAGS  . 

Compare  them  with  Figures  0.4,  D.ll,  D.13,  and  D.16. 

The  magnitudes  of  the  FWT  coefficients  at  each  In]  are  not 
the  same  in  the  shifted  TAGS.  Compare  TOA  TAG  S  and  TOA  TAG  S 
SHF  1.  However,  notice  that  shifted  TOA  TAG  12  <T0A  TAG  12 
SHF  1)  has  the  same  magnitude  at  each  WALSH  index  Cnl  as  does 
the  unshifted  TOA  TAG,  but  the  components  from  n  *  32  to  n  = 

63  are  opposite  in  sign.  The  even  symmetry  of  the  even 
numbered  FWT  about,  the  midpoint  has  been  changed  to  odd 
symmetry  about  midpoint. 

This  is  a  very  neat  and  clean  distinction,  and  further 
emphasizes  the  distinction  between  even  and  odd  numbered  TOA 
TAGS.  It  is  further  evidence  that  at  least  some  general 
features  exist  which  could  be  made  to  work  in  a  computer 
algorithm  to  identify  individual  PRI's. 

FWT's  of  odd  numbered  intervals  have  no  general 
distinctive  change  between  shifted  plots  and  standard  plots. 

In  a  few  shifted  odd  numbered  plots,  there  are  some  particular 
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Figure  4.2.  FWT  j  TOA  TAG  5  SHIFT  I 
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magnitudes  o-f  index  Cn]  that  remain  the  same  or  increase 
proportionately,  but  others  do  not.  If  any  real  relationship 
between  the  shifted  and  unshifted  odd  numbered  FWT  plots 
exist,  then  it  will  have  to  be  revealed  with  a  much  more 
detailed  study. 

Quite  simply,  odd  numbered  shifted  FWT  plots  show  no  major 
distinctive  features  that  are  readily  apparent.  Even  numbered 
shifted  FWT  plots  follow  the  sign  change  rule  for  components 
at  n  »  32  to  63,  but  the  magnitudes  of  the  components 
of  both  plots  are  the  same. 

Although  the  odd  numbered  PSD  plots  change  in  the 
amplitude  of  the  sequency  components  from  unshifted  to 
shifted,  the  even  numbered  PSD^s  ar  •  the  same  in  both  shifted 
and  unshifted  plots.  Indeed,  they  are  the  same  in  both 
magnitude  and  location  of  sequency  components.  Compare 
Figures  4.6  to  4.9  on  the  following  pages  with  Figures  D.26, 
D.33,  D.3S,  and  D.38  in  Appendix  D. 

Thus,  the  even  numbered  PSD's  of  this  DATA  line 
representation  join  the  time  invariant  Group  Spectrums 
in  being  immune  to  the  effects  of  a  circular  time  shift. 
Compare  the  shifted  Group  Spectrums  Figures  4.18  to  4.13,  with 
the  unshifted  Group  Spectrums  in  Appendix  E,  Figures  E.l 
through  E.4. 

Again  the  point  is  made  that  these  general  features  could 
be  recognized  in  a  PRI  recognition  algorithm.  At  the  very 
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lea&t,  it  Mould  eliminate  the  odd  or  even  choices  o-f  time 
intervals  between  numbers  of  the  single  PRI  representation. 

D,  INTERLEAVED  TOA  STRINGS  AND  THE  FWT 

The  receipt  of  radar  pulses  from  several  different 
emitters  results  in  individual  pulse  trains  that  are 
interleaved  into  a  data  stream  from  the  receiver.  Separation 
or  deinterleaving  of  these  pulse  trains  from  the  data  stream 
must  be  done  to  identify  the  individual  emitters.  This 
section  examines  the  plots  of  the  FMT  and  PSD  coefficients  of 
a  DATA  line  representation  of  interleaved  pulse  trains 
composed  of  the  TOA  TAG  representation  in  Section  B.  of  this 
chapter. 

It  is  desired  to  answer  the  following  questions: 

1.  Are  the  components  of  the  FWT  and  PSD  of  an 
interleaved  representation  identifiable  as  belonging  to  a 
particular  PRI? 

2.  Is  there  a  relationship  between  the  magnitude  or 
sequency  of  the  components  of  an  interleaved  plot  with  that 
of  plots  of  single  PRI  trains? 

3.  Are  interleaved  plots  made  by  simply  adding  single  PRI 
pi ots? 

Begin  by  considering  interleaved  TOA  TAGS  of  Table  4.2. 
The  DATA  line  representation  for  TOA  TAG  5  with  N  «  64  is 
DATA  IjdiOfOyOiliBiOfOySilfOfSySfOyljOiBiBiSflyOiGfBfGilydi 
0, 1,0, 0,0,0, 1,0, 0,0  <4.3> 

and  for  TOA  TAG  8, 

DATA  1,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 1,0,0, 
0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0,0, 
0,0, 1,9, 0,0, 0,0, 0,0  (4.4) 
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Now,  i n ter 1  save  the  two, 

i,e,e,e,8,i,e,e,i,e,i,e,e,e,e,i,i,0,8,e,i,e,0,8,i,t,0, 

8, 0,0, 1,8, 1,8, 8, 1,0, 0,8, 0,1, 8, 8, 8, 8, 1,8, 8, 1,0, 1,0, 8, 8, 

8, 1,  1,8, 8, 9, 1,8, 8,8  (4,5) 

Notice  that  the  -i^irst  pul  so  o-f  each  TAG  is  coincident  and 
both  are  represented  by  the  -first  pulse  o-f  the  interleaved 
DATA  line. 

The  FVrr  o-f  TAGS  3  and  8  are  shown  in  Figures  4.14  and 
D.7.  TOA  TAG  3  components  are  emphasized  by  setting  the 
WAL(0,t}  components  equal  to  8  before  plotting.  The  plotting 
routine  then  normalizes  to  a  smaller  magnitude  component  and 
increases  the  relative  magnitude  of  the  other  components. 

Examine  the  interleaved  FWT,  Figu'^e  4.  IS.  Note  sizable 
components  at  9,  22,  25,  31,  32,  37,  45,  48,  58,  60,  and  63. 
Now  examine  TOA  TAG  3  with  WAL<0,t>  equal  to  8.  Note 
emphasized  components  at  25,  37,  44,  45,  50,  51,  and  52. 

There^s  some  commonality  between  the  plots.  Components  at 
25,  37,  45,  50,  and  52  are  obviously  associated  with  TOA  TAG 
5.  So  where  is  TOA  TAG  8?  A  closer  1 ooK  reveals  its 
influence  on  the  interleaved  FWT.  Note  identical  relative 
magnitude  pairs  at  n  ~  IS  and  16  and  n  »  31  and  32.  These  are 
generated  by  the  coefficient  pairs  of  these  Cnl  from  TOA  TAG 
8.  Note  also  that  n  *  48  of  the  interleaved  FVfT  is  highe.'  in 
relative  magnitude  because  there  is  a  small  n  »  48  component 
of  TOA  TAG  5. 

How  would  one  use  this  interleaved  FVir  to  deinterleave 
pulse  trains?  A  direct  computation  method  using  a  correlation 
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algorithm  is  considered  in  the  next  section.  The  author 
suggests  a  ‘matching"  algorithm  could  be  written  that  would 
scan  the  interleaved  FWT,  and  note  relatively  high  magnitudes 
at  certain  Cn]  above  some  threshold.  These  tnl  could  be 
matched  to  a  data  bank  o-f  important  Cnl  components  o-f  single 
PRI  FWT  plots,  along  with  any  special  characteristics  o-f  the 
interleaved  FWT  < identical  components  at  successive  Cnl,  for 
example).  Promising  matches  between  the  interleaved 
parameters  and  single  PRI  parameters  would  result  in  a 
decision  or  choice  of  which  single  PRI  trains  make  up  the 
interleaved  pulse  train. 

Scrutiny  of  the  PSD  plots  generates  similar  conclusions. 
Figure  4.16  shows  the  interleaved  PSD  plot.  Components  of 
relative  magnitude  of  0.3  or  greater  are  8,  13,  16,  19,  22, 

24,  25,  26,  and  32. 

Now  note  the  PSD's  of  TOA  TAG  5  and  8  (Figures  4.17  and 
D.29w  TOA  TAG  5  contributes  components  at  sequencies  of  13, 
22,  25,  and  26.  TOA  TAG  8  contributes  components  at  8,  16,  24, 
and  32. 

Before  moving  on  to  another  interleaved  case,  a  discussion 
of  the  magnitudes  of  the  interleaved  components  (both  FWT  and 
PSD)  is  in  order.  A  definite  relationship  is  not  apparent 
between  components  of  individual  plots  and  interleaved  plots. 

The  author  can  say  that  one's  intuition  about  components 
of  interleaved  plots  generally  holds  true.  For  instance,  a 
large  positive  FWT  component  of  one  TAG  and  a  negative  FWT 
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Figur*.  4.16.  PSD  :  TOA  TAG  3  FWD  8  INTERLEAVED. 
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component  o-f  another  TAG  at  the  same  tn3  will  usually  generate 
a  smaller  component  in  a  FWT  component  o-f  the  interleaved 
TAG's^  The  result  is  not  a  simple  addition,  but  the  actual 
magnitude  o-f  the  interleaved  component  is  usually  smaller  than 
the  added  values  o-f  the  single  components.  No  fast  and  easy 
rule  appi ies. 

In  a  matching  algorithm,  the  magnitude  would  have  only  a 
comparative  bene-fit.  The  actual  promise  of  the  algorithm  lies 
in  matching  component  locations  <either  Cnl  or  sequency)  to  a 
data  bank  o-F  single  PRI  component  locations,  with  magnitude 
and  other  special  characteristics  as  a  decision  making  aid. 

One  more  interleaved  case  is  examined.  This  one  is 
di-f-ferent  because  there  are  no  coincident  pulses.  The  shift 
in  TQA  TAG  3  causes  it  to  interleave  perfectly  with  TOA  TAG 
19. 

TGA  TAG  5  SHF  1  <Figure  4,2  and  4.4)  DATA  line: 

DATA  0, 1,0, 0,0,0, 1,0, 0,0,0 , 1,0, 0,0,0, 1,0, 0,0,0 , 1,0, 0,0,0, 1,0, 

0,0,0, 1,0, 0,0,0, 1,0, 0,0, 0,1, 0,0, 0,0, 1,0, 0,0,0, 1,0, 0,0, 

0, 1,0, 0,0,0, 1,0,0  (4.5) 

TQA  TAG  19  (Figure  4.18  and  4.19)  DATA  line; 

DATA  1,0, 0,0, 0,0, 0,0, 0,0, 8, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,0, 
0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0, 1,0, 0,0, 0,0,0  (4.6) 

And  the  interleaved  DATA  line; 

DAT  A  1,1, 0,0, 0,0, 1,0, 0,0, 0,1, 0,0, 0,0, 1,0, 0,1, 0,1, 0,0, 0,0, 1,0, 

0,0, 0,1, 0,0, 0,0, 1,0, 1,0, 0,1, 0,0, 0,0, 1,0, 0,0, 0,1, 0,0,0, 

0, 1, 1,0, 0,0, 1,0,0  (4.5) 

Figure  4.20  and  4.21  show  the  FWT  and  PSD  of  the 

interleaved  TAG'S.  Interleaved  components  of  both  single  PRI 
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TAGS  are  present  v.  and  thi?  same  process  as  be-fore  is  used  to 

( 

Identi-fy  them.  A  he*4  discovery  of  this  examination  reveals 
that  some, feature  PSD  components  of  TOA  TAG  5  SHF  I  are  passed 

directly  to  the  i r t er 1 oaved  PSD.  Note  the  magnitude  of 

/ 

sequencies  'll,  16,  and  interleaved  PSD  is  the  same 

as  the  PSD  of  TOA  TAG,  Z  SHF  Ij  and  that  the  magnitudes  of  TOA 
TAG  17  PSD  wt  thesh  scqueticies  is  xero.  This  is  the  first 
evidence  of  two  single  train  representations  generating  an 
inter  1  eav;>d  components  magnitude  by  merely  adding  the  simple 
PHI  components.  A  more  c]v.<snr  study  could  be  done  to  see  if 
it  is  merely  coincidence.  The  real  question  is,"  Did  the  fact 
that  there  are  no  coiricident  pulses  in  the  interleaved 
represen tation  allow  me  to  add  certain  single  PRI  components 
to  form  an  interleaved  component?" 

A  partial  answer  to  the  question  probably  is  evident  since 
only  certain  components  of  the  two  single  PRI  PSD  plots  have 
this  feature  of  ‘dimple  addition. 

E.  coRnrLATinri  pnd  r;t£  fSD  coefficients 

Throughout  ^his  chapter  the  emphasis  has  been  placed  on 
the  1  oration  and  relative  magnitudes  of  components  located  in 
the  Malsh  index  domain  or  the  sequency  domain.  Most  of  the 
plots  show  the  magnitude  of  the  components  relative  to  the 
maximum  component  magnitude  of  that  particular  plot.  It  was 
stated  several  times  that  the  actual  magnitudes  of  the 
individual  PRI  <TOA  TAG)  components  didn't  seem  to  have  a 
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recognizable  relationship  with  the  magnitudes  of  the 
components  of  the  interleaved  plots. 

It  was  thought  that  the  magnitudes  might  be  related,  and  a 
correlation  of  the  interleaved  component  magnitudes  with  those 
of  the  individual  PRl  magnitudes  might  show  their  presence. 

The  PSD  coefficients  were  chosen  because  they  are  defined  from 
the  FWT  coefficients,  and  exhibit  the  same  characteristics  and 
information  about  the  transform  in  terms  of  symmetry  and 
relative  magnitude  components. 

Without  any  mathematical  fanfare,  the  correlation 
coefficient  p  is  presented. 

P  =  E  <  <X-XO/<J5^  •  (Y-Y^)/dY  >  <4.6> 

where 

X,Y  are  random  variables 

X^,  Y''  are  the  mean  of  the  variables 

variances  of  the  variables 

The  correlation  coefficient  expresses  the  degree  to  which 
two  random  variables  are  correlated  without  regard  to  the 
magnitude  of  either  one.  C30:  p.83]  It  is  a  normalized 
quantity  and  will  be  a  value  between  -'I  and  +1. 

In  this  application,  X  are  the  PSD  coefficients  of  the 
interleaved  situation  and  Y  are  the  PSD  coefficients  of  the 
individual  TAG  situation. 

A  correlation  BASIC  language  routine  was  written 
specifically  for  the  32  PSD  components  of  the  interleaved 
and  individual  TOA  TAGS.  This  program  was  adapted  from  a 
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correlation  routine  in  the  IMSL  library  contained  in  the  Naval 
Postgraduate  School  IBM  3833  computer.  It  was  not  adapted  for 
general  application  but  specifically  for  this  situation.  For 
this  reason  it  is  not  provided  in  Appendix  A. 

The  correlation  routine  merely  reads  in  the  component 
values  of  th<|>  individual  PRI  TAG  representations  and  the 
interleaved  representation  generated  by  the  FWT/PSD/Group 
Spectrum  program  provided  in  Appendix  A.  It  correlated  the 
two,  based  upon  Equation  (4.d> ,  and  outputed  a  single  number, 
the  correlation  coefficient. 

The  PSD  coefficients  of  the  interleaved  situation  of  TOA 
TAG  S  and  TOA  TAG  19  were  used  < Figure  E.S  in  Appendix  £> . 
This  interleaved  situation  is  different  from  that  of  page  148 
because  TOA  TAG  5  wasn't  shifted  before  interleaving,  causing 
the  first  pulse  of  TOA  TAG  5  and  19  to  be  coincident.  Nith 
these  coefficients,  the  values  of  the  individual  TOA  TAG  ^ 
components  from  3  to  64  were  then  correlated  with  the 
interleaved  values  using  the  correlation  routine. 

Table  4.2  gives  the  TOA  TAG  number,  and  the  correlation 
coefficient  generated. 

Examination  of  the  correlation  coefficients  shows  strong 
correlation  with  TOA  TAG  5.  Great!  It  is  desired  to  have  a 
high  correlation  value  with  TOA  TAG’S  3  and  19,  since  the 
interleaved  representation  was  composed  of  these  two  TAG 
representations.  The  TOA  TAG  5  correlation  coefficients  fits 
the  need  nicely. 
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Table  4.2.  Correlation  o-f  Interleaved  PSD  coe-f-f icients  uiith 

individual  TOA  TAGS 


TOA  TAG  number 


Correlation  Coe-f-f icient 


3 

0.665 

4 

0,234 

5 

0.963 

6 

0.582 

7 

0.305 

8 

0.030 

9 

0.166 

10 

0.709 

11 

0.2^9 

12 

0. 183 

13 

0.274 

14 

0. 191 

15 

0.365 

16 

-0.065 

19 

0.140 

22 

0.330 

23 

-0.028 

32 

-0.055 

40 

0 . 106 

51 

0.059 

64 

-0.722 

On  the  down  side,  the  interleaved  PSD  coe-f-f  icients  do  not 
show  a  high  value  o-f  correlation  with  TOA  TAG  IP.  The  initial 
success  o-f  TOA  TAG  5  certainly  indicates  that  this  is  an  area 
for  further  investigation,  however. 

The  author  suggests  that  the  correlation  procedure  could 
be  used  in  an  effective  deinterleaving  algorithm  if  it  was 
applied  to  the  sequential  method  of  extracting  PRI^s.  Once  the 
initial  correlation  routine  is  carried  out  and  a  high  value  of 
correlation  coefficient  is  identified  and  associated  with  an 
individual  TOA  TAG  string,  then  the  algorithm  could  return  to 
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the  interleaved  DATA  stream  representation  and  subtract  <set 
equal  to  zero)  those  TOA  TA6‘'s  associated  with  the  determined 
PRI . 

Now  repeat  the  entire  process.  TaKe  the  FWT  and  PSD  of 
the  remaining  TOA  TAG  representation,  correlate  the  PSD 
coefficients  with  the  data  bank  of  individual  TOA  TAG  PSD 
coefficients,  and  see  if  a  new  TOA  TAG  has  been  identified 
with  another  high  value  of  correlation  coefficient.  If  it 
has,  then  strip  these  TOA  TAG'^s  from  the  original  string. 

The  process  could  be  repeated  until  calculation  of  the  PSD 
yields  no  further  useful  coefficients,  or  until  no  correlation 
IS  found.  The  stripped  TAG'S  represent  the  individual  pulse 
trains  that  make  up  the  interleaved  pulse  string. 

This  suggested  process  is  similar  to  the  way  a 
preprocessor  sorts  the  incoming  data  stream.  After 
identifying  particular  pulses  in  the  stream,  the  preprocessor 
strips  these  pulses  from  the  stream,  simplifying  the  follow-on 
process,  which  continues  to  identify  <deinter 1 eave)  individual 
pulse  trains  of  the  remaining  stream.  Pach  identified  train 
is  then  removed  from  the  incoming  data  from  the  receiver. 

A  problem  that  plagues  current  deinterleaving  techniques, 
that  of  coincident  or  missing  pulses,  would  probably  plague 
this  correlation  deinterleaving  procedure. 

The  TOA  TAG  representation  does  not  indicate  coincident 
pulses  separately.  Coincident  pulses  are  represented  by  a 
single  TOA  TAG.  Only  the  receiver  knows  which  one  it  detected 
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and  rspresantad  as  a  received  pulse.  Therefore,  the  stripping 
of  this  TAG  generated  when  two  pulses  had  arrived  at  the 
receiver  at  the  same  time  inadvertently  strips  a  TOA  TAG  from 
the  remaining  pulse  representation. 

If  the  PSD  coefficients  of  the  remaining  string  with  the 
inadvertently  stripped  TAG  were  now  calculated,  one  would  find 
they  are  different  than  those  calculated  from  the  remaining 
string  that  included  the  inadvertently  stripped  TAG. 

These  PSD  coefficients,  those  calculated  from  the  missing 
TAG  string,  cannot  correlate  higher  with  the  individual  TOA 
TAG  strings  <used  by  the  correlation  routine)  than  the 
coefficients  of  the  complete  remaining  string. 

The  individual  TAG  PSD  coefficients  were  calculated  from 
TOA  TAG  strings  with  no  missing  pulses.  So,  coefficients  from 
strings  with  missing  pulses  will  not  correlate  as  well  as 
coefficients  calculated  from  complete  strings, 

A  possible  solution  could  be  the  use  of  a  large  number  of 
TOA  TAG'S  in  the  interleaved  representation,  say  1024  or  2048. 
In  other  words,  expand  the  "window*  or  aperture  on  the  data 
stream.  This  would  allow  the  interleaved  representation  to  be 
more  realistic  and  comparable  to  the  receiver's  own  TOA 
resolution.  More  coefficients  might  help  the  correlation 
procedure  and  lessen  the  effect  of  missing  and  coincident 
pu 1 ses . 
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V.  CONCLUSIONS  AND  RECOhMENDAT I ONS 


B«-fore  discussing  the  conclusions  of  the  research  effort 
and  maKing  recommendations  for  follow-'on  work,  a  review  of  the 
investigation  steps  that  were  taken  by  the  author  and 
translated  into  the  substance  of  this  thesis  will  be 
beneficial . 

For  background,  the  importance  of  the  ESM  effort  and  its 
place  in  Electronic  Warfare  was  studied.  A  typical  ESN 
system  was  examined,  and  the  problem  of  deinterleaving  the 
pulse  stream  output  from  an  ESN  receiver  was  considered,  along 
with  current  techniques  for  deinterleaving  already  in  use. 

Having  no  previous  idea  of  even  the  existence  of  the  Walsh 
functions  set  the  stage  for  the  study  of  these  unique  and 
interesting  functions.  As  the  study  proceeded  through  the 
definition  of  the  Walsh  Transform,  the  constant  comparisons  of 
the  Walsh  Transform  properties  with  the  familiar  properties  of 
the  Fourier  Transform  brought  review  of  this  area. 

Reference  23  provided  the  FWT  program  that  became  the  too] 
of  the  investigation.  Adaptation  of  the  program  to  the  BASIC 
language  and  a  gathering  of  evidence  that  supported  belief  in 
its  output  took  a  considerable  time.  Searches  of  the 
literature  for  examples  of  FWT''s  of  familiar  waveshapes 
finally  provided  the  proof  of  the  program. 

With  the  guidance  of  the  thesis  advisor,  the 
representation  of  the  receiver  data  stream  with  time  of 


Arrival  tags  was  accepted,  and  the  FKT  program  was  applied  to 
many  individual  tag  strings.  The  purpose  was  to  determine  i-f 
a  TOA  TAG  string  representing  a  particular  PRI  pulse  train 
could  be  recognized  by  certain  -features  o-f  the  Mai  sh  Trans-form 
or  the  Power  Spectral  Density.  The  FUTT,  PSD,  and  Group 
Spectrum  o-f  interleaved  TOA  TAGS  were  then  closely  examined 
for  PRI  recognition  features.  Finally,  a  correlation 
algorithm  was  used  to  determine  the  degree  of  correlation 
between  interleaved  PSD  coeffcients  and  single  pulse  train  PSD 
coefficients.  An  examination  of  the  correlation  coefficients 
was  done  to  see  if  a  PRI  could  be  determined  by  noting  the 
magnitude  of  the  coefficient. 

A.  CONCLUSIONS 

One  result  of  the  effort  was  the  demonstration  of  the 
ease  of  computation  of  the  Walsh  Transform  and  Power  Spectral 
Density  of  a  number  series.  The  algorithm  is  fast  and 
efficient,  and  easily  adapted  to  a  SASIC  language 
implementation.  This  allows  the  computation  of  the  transform 
and  calculation  of  the  PSD  coeffcients  on  a  microcomputer, 
easily  used  and  readily  available  in  the  worK  space. 

When  applied  to  a  TOA  TAG  representation  of  a  single  PRI 
pulse  train,  the  transform  snows  features  that  can  be  used  to 
recognize  this  pulse  train  when  compared  to  the  transforms  of 
other  PRI  pulse  trains.  Symmetry  of  the  coefficient 
distribution,  and  the  relative  magnitude  of  certain  Walsh 
index  or  sequency  components  are  the  features  that  allow  this 
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identification 


It  was  shown  in  Chapter  4  that  each  single  PRI  pulse 
train's  Walsh  transform  has  unique  and  distinctive  components 
that  distinguish  it  from  other  transforms  of  single  PRI 
representations.  No  single  feature,  however,  such  as  a 
relative  maximum  sequency  component,  can  be  connected  with  the 
PRI  of  the  pulse  train.  Each  particular  PRI  transform  simply 
has  individual  features  th.*t  identify  it,  and  only  certain  ’ 
characteristics,  such  as  symmetry,  are  common  between  groups 
of  individual  pulse  train  transforms. 

The  Walsh  transform  and  PSD's  of  interleaved  TOA  TAG 
representations  have  features  that  can  be  identified  as 
belonging  to  the  transforms  or  PSD's  of  component  individual 
pulse  trains.  However,  these  features  are  not  readily 
apparent  and  one  could  not  identify  or  associate  them  with  a 
particular  single  PRI  transform  without  prior  Knowledge  of 
single  PRI  transform  features.  LiKe  the  transforms  of  single 
PRI  pulse  trains,  the  interleaved  transform  or  PSD  does  not 
have  any  features  that  can  identify  component  PRI's  from  the 
transform  or  PSD  alone. 

This  lack  of  component  PRI  features  in  the  transform  or 
PSD  effectively  prevents  its  use  in  a  practical  deinte»'l  eaving 
algorithm.  At  the  outset  of  the  research  effort,  it  was  hoped 
that  the  transform  of  a  representation  of  the  pulse  stream 
from  a  ESM  receiver  would  identify  the  individual  PRI's  of  the 
component  pulse  trains.  Neither  the  transform  or  the  PSD  does 
Inis,  al though  there  are  common  features  between  the 
interleaved  transform/PSD  and  component  transform/PSD. 


The  application  o-f  a  correlation  routine  between  the 
the  interleaved  PSD  coefficient's  and  the  individual  TOA 
TAG  PSD  coefficients  was  seen  to  be  effective  in  identifying 
one  of  the  cotnponent  PRI''s.  It  did  not  identify  all 
component  PRI's,  and  any  algorithm  that  uses  this 
approach  will  have  to  go  through  additional  steps  to  do  so. 
These  required  additional  steps  quickly  negate  the  speed  of 
computation  advantage  that  the  Hal sh  transform  exhibits;  and 
the  approach  becomes  quite  similar  to  current  sequential  PRl 
extraction  techniques. 

In  summary,  the  transform  of  the  interleaved  TOA  TAG 
representation  failed  to  produce  a  simple  PRI  recognition 
feature.  Recognition  features  are  present,  however,  in  the 
interleaved  transform  and  PSD  that  would  allow  the 
identification  of  component  pulse  trains  in  an  interleaved 
representation  if  a  data  bank  of  individual  PRI  transform/PSD 
features  was  available. 

B .  RECOftlENDAT  I  QMS 

It  is  recommended  that  an  effort  at  writing  a  "matching 
feature*  deinterleaving  program  be  made.  The  program  would 
use  a  data  bank  of  parameters  that  describe  the  interesting 
properties  and  features  of  the  Walsh  transform  and/or  PSD  of 
TOA  TAG  representation  of  a  single  PRI  pulse  train. 

The  data  bank  would  be  the  measure  of  attaining  a 
successful  deinterleaving  program.  Properties  such  as  even 
symmetry,  double  components  of  power  of  two  interval 


representations,  and  distinctive  index  or  sequency  components 
o-f  the  odd  numbered  TOA  TAG  representations  would  be 
inci uded. 

The  algorithm  would  scan  the  FWT/PSD  o-f  the  interleaved 
TOA  TAG  representation  for  recognizable  features,  then  match 
these  f<».Ttures  to  the  data  bank.  Appropriate  tolerances  could 
be  included  to  be  used  in  decisions  about  whether  a  single  PRI 
is  a  component  of  the  interleaved  train. 

Additional  study  could  possibly  yield  the  optimal  features 
to  include  in  the  bank.  For  certain,  a  closer  look  at  single 
PRI  representations  must  be  done  to  determine  which  components 
are  shared  among  them.  The  author  firmly  believes  that  there 
will  always  be  a  FWT/PSD  feature  of  a  single  PRI 
representation,  whether  a  particular  index  or  sequency 
component,  or  a  symmetrical  characteristic,  that  can  be  used 
to  distinguish  it  from  other  single  PRI  representations. 

Closer  study  of  whether  there  is  a  relationship  between 
interleaved  and  individual  index  or  sequency  magnitudes  should 
be  completed.  It  is  a  feature  that  could  possibly  be  used  in 
a  matching  algorithm,  and  it  certainly  can  be  used  in  direct 
compuation  methods,  such  as  a  least  mean  squares  or 
correlation  routine  that  might  generate  an  individual  PRI 
indication. 

The  quite  tentative  but  interesting  result  found  in  the 
last  section  of  Chapter  4  should  be  explored  completely.  All 
coefficients  should  be  involved,  including  the  Group  Spectrum 
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coe-f-f icients.  In  this  thesis,  the  Group  Spectrum  coefficients 
have  mainly  been  used  to  show  the  time  invariance  nature  of 
their  character,  but  they  might  offer  some  interesting  direct 
computation  features. 


APPENDIX  A 


FAST  WALSH  TRANSFORM  PROGRAMS 

The  -following  pages  contains  listings  o-f  the 
Walsh  Trans-form  programs  used  in  this  thesis.  The  -first 
program  is  written  for  an  IBM  Personal  Computer,  and  the 
second  was  written  for  the  HP-85. 

The  FWT  subroutine  was  adapted  from  a  FORTRAN  program  in 
"Walsh  Functions  and  their  Applications,”  by  K.6.  Beauchamp. 

The  PSD  and  Group  Spectrum  subroutines  were  written  by  th 
author . 

Both  programs  were  run  on  the  respective  computer  BASIC 
interpreter  < uncompi 1 ed> .  A  compiled  version  would  probably 
run  faster,  although  a  128  number  FWT  ran  on  the  IBM  PC  in 
about  4  seconds  using  only  the  interpreter. 
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10  '  This  program  computes  the  Fast  Walsh  Tranform  ,  the 
Power  Spectral  Density,  and  the  Group  Spectrum 
20  '  Coe-f-f icients  o-f  a  number  series. 

30  '  The  FWT  subroutine  is  a  BASIC  adaptation  o-f  a  FORTRAN 
program  listed  40  ^  in  "  Walsh  Functions  and  Their 
Applications,*  by  K.6.  Beauchamp. 

45  '  The  PSD  and  Group  Spectrum  subroutines  were  written  by 
the  author. 

50  '  This  program  runs  on  a  IBM  Personal  Computer.  The 
printer  used  was  a  NEC  8023. 

60  DIM  X( 100) ,Y<64) ,P<64) ,Q< ll> ,F(64)  'dimension  arrays 
65  LPRIMT  CHR»<27)  j*L- ;"008"  'set  le-ft  margin 

70  N  s»  64  '  number  o-f  input  numbers 

80  FOR  1=1  TO  N  '  read  in  data  ......... 

90  'READ  X<I)  '-from  the  DATA  line,  or........... 

100  X<I)  =  SIN( 1*I/N*2*3. 1416)  'calculate  your  own 
110  NEXT  'I 

120  GOSUB  310  'calculate  FWT  o-f  input  series 

130  FOR  B  =  1  TO  N  'divide  a1  1  coe-ff s  by  N  and  downsh-ft  1 
140  F<B-1)  »  X<B)/N  '  this  array  contains  the  FWT  coe-ffs 
150  NEXT  'B 

169  GOSUB  600  'calculate  power  spectrum  coe-f -f i c i en ts 

170  GOSUB  720  '  calculate  group  spectrum 

180  GOSUB  840  'output  the  FWT  coe-ffs 

190  GOSUB  950  'output  the  PSD  coeffs 
200  GOSUB  1047  'output  the  group  coeffs 
230  END 
235  ' 

240  '  SUBROUTINE  FHT<N,X,Y) 

250  ' 

260  '  This  routine  performs  a  FWT  of  an  input  series  in  array 
X.  The  array  Y  is  used  for  working  space. 

270  '  The  dimensions  of  X  and  Y  must  be  a  power  of  2. 

280  '  The  results  of  the  FWT  are  in  sequency  order,  positive 
290  '  phasing,  and  in  array  X.  The  subroutine  uses  a  Hadamard 
transform. 

308  ' 

310  N2  =  N/2 

320  M  ■  6  '  M  equals  the  base  2  logarithm  of  N 

330  FOR  L  -  1  TO  M 


340 

NY  =  0 

350 

NZ  »  2 

*<L 

-I) 

360 

NZI  = 

2  X 

NZ 

370 

NZN  » 

N/NZI 

380 

FOR  I 

»  1 

TO 

NZN 

390 

NX  » 

NY 

1 

400 

NY  = 

NY 

-t- 

NZ 

410 

JS  = 

<I 

-1) 

X  NZI 

420 

JD  = 

JS 

NZI  ♦ 

430 

FOR 

J  = 

NX 

TO  NY 

440 

JS 

— 

JS 

-►  1 
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450  J2  “  J  +  N2 

460  Y<JS>  =«  X(J)  +  X<J2> 

470  JD  *  JD  -  1 

480  Y<JD>  =  X(J>-  X<J2) 

490  NEXT  "J 

500  NEXT  ' I 
510  FOR  B  =  1  TO  N 
520  X<B)  =  Y<B) 

530  NEXT  '9 
540  NEXT  'L 
550  RETURN 

555  '  SUBROUTINE  PSD<N) 

560  ' 

570  ' 

580  '  This  routine  calculates  the  Nalsh  Power  Spectral 
Density  Coefficients. 

590  ' 

600  P<0>  *  <F<0))^2  'first  PSD  coefficient 

610  K=1 

620  FOR  I  »  1  TO  <N/2)-l  'coefficients  1  to  N/2-1 
630  P<l)a«  <F<K))'‘2  +  <F<K+1))'‘2 

640  K  *  K  +  2 

650  NEXT  'I 

660  P<N/2)*<F<N-i))'‘2  'last  PSD  coefficient 

670  RETURN 
<^80  ^ 

690  '  SUBROUTINE  GROUP  SPECTRUMCM) 

700  ' 

719  '  This  subroutine  calculates  the  group  spectrum 

coefficients. 

720  D»1 

730  Q<0)  -  P<0)^2 
740  FOR  8  ■  0  TO  M-l 

750  FOR  C  =»<2'‘B)  TO  N/2-1  STEP  (2X2'‘B> 

760  G<0)  »  6<D>  +  P<C)*2 

770  NEXT  'C 

780  D  *  0+1 

790  NEXT  '8 

800  S<M)  -  P<N/2)*2 

810  RETURN 

820  ' 

830  '  This  subroutine  outputs  the  Walsh  Transform 
Coefficients. 

840  ' 

841  FOR  C  »  1  TO  5:  LPRINT:  NEXT  '  1  inch  top  margin 

842  LPRINT*  FWT  COEFFICIENTS  FOR  1  CYCLE  SINE 

FUNCTION* 

843  LPRINT;  LPRINr  'skip  two  lines 

850  LPRINT  'Walsh  Function'  TAB<20)  'Coeff."  TAB< 35> *Wal sh 
Function'  TAB<55>  'Coeff,' 
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860  LPRIhTT  STRINe*<  15,45)  TA8<20)  STRING*<  6 , 45)  TAB(35) 
STRING*< 15,45)  TAB<55)  STRING*<6,45) 

878  FOR  8=1  TO  32  "output  the  coe-f-fs 

880  LPRINT  "mLC  jB-l;“  ,t)  •  TA8(  16)  :  LPRINT  USING 

5F<8-1)  ;  !  LPRINT"  ■  TA8<35) 

890  LPRINT  •WAL<" ;B+31i" ,t) ■  TAB<52)  :  LPRINT  USING 

;F(B4-31) 

900  NEXT  '8 

910  FOR  C  =  1  TO  16  ;  LPRINT:  UOCT  "skip  16  lines 

911  LPRItTT"  xx"  "print  pg.  no 

912  LPRINT  CHR»<12)  "-form  -feed 
920  RETURN 

930  " 

940  "  Output  the  PSD  coefficients. 

950  " 

951  FOR  C  =  1  TO  5;  LPRINT:  NEXT 

970  LPRINT  ■  PSD  COEFFICIENTS  FOR  A  1  CYCLE  SINE 

FUNCTION" 

980  LPRINT  :  LPRINT 
990  FOR  8  =  0  TO  N/2 

1000  LPRINT  "P<";  8  ;")="  LPRINT  USING  "MM .  ;  P<  8) 

1020  NEXT  "B 

1025  LPRINT:  LPRINT 

1040  RETURN 

1044  " 

1045  '  Output  the  group  coefficients 

1046  " 

1047  LPRiNT"  GROUP  SPECTRUM  COEFFICIENTS  OF  A  I  CYCLE  SINE 

FUNCTION" 

1048  LPRINT:  LPRINT 
1050  FOR  B  »  0  TO  M 

1060  LPRINT  •G<"jBj")»  "  }  :LPRINT  USING  ;G<B) 

1070  NEXT  "B 

1071  FOR  C  »  1  TO  7:  LPRINT:  NEXT 

1072  LPRINT  •  XX" 

1073  LPRINT  CHR»< 12)  'form  feed 
1080  RETURN 
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This  program  computes  the  Fast  Mai  sh  Tran-form  ,  the 
Power  Spectral  Density,  and  the  Group  Spectrum 

29  *  Coefficients  of  &  number  series. 

30  i  The  FWT  subroutine  is  a  BASIC  adaptation  of  a  FORTRAN 

program  listed  40  !  in  *  Nalsh  Functions  and  Their 
Applications,”  by  K.G.  Beauchamp. 

45  !  The  PSD  and  Group  Spectrum  subroutines  were  written  by 
the  author. 

50  This  program  runs  on  HP  85  Computer,  and  writes  the 

coefficients  to  data  files  named  ”  Fl^”  ,  ■  PSD*  ,  and  "GRP". 

69  DIM  X(  10?>  ,Y<64>  ,P<64)  ,G<  11)  ,F<64)  '.dimension  arrays 

65  PRINTER  IS  2  !  select  printer  for  output  of  print 

70  N  a  64  <  number  of  input  numbers 

89  FOR  I  s  1  TO  N  •  read  in  data  . 

90  !REAO  X<I)  ‘from  the  DATA  line,  or . . . 

100  X<I)  =  SIN< IXI^NXaxa. 1416)  ‘calculate  your  own 
110  NEXT  I 

120  GOSUB  310  ‘calculate  FWT  of  input  series 

130  FOR  B  =  1  TO  N  ! divide  all  coeffs  by  N 

140  X<B)  »  X<B)/N  !  this  array  contains  the  FWT  coeffs 

150  NEXT  B 

155  BEEP  e  CLEAR 

156  DISP  "FWT  CAL  OUER" 

160  GOSUB  696  ! calculate  power  spectrum  coefficients 

165  DISP  "PSD  CAL  OVER” 

170  GOSUB  720  !  calculate  group  spectrum 

175  DISP  "GRP  COEFF  CAL  OVER- 
188  Q0SU6  840  ‘output  the  FWT  coeffs 

190  GOSUB  950  ! output  the  PSD  coeffs 

200  GOSUB  1047  ! output  the  group  coeffs 

205  CREATE  -FWT”, 30  !  creates  file 

206  ASSIGN  Ml  TO  ”FWT'  !  opens  file 

207  FOR  8  =  1  TO  64 

208  PRINTM  1}  X<0) 

209  NEXT  0 

210  ASSIGNM  1  TO  X  ! Cl oses  file 

211  CREATE  ”PSD”,  26 

212  ASSIGN  Ml  TO  "PSD' 

213  FOR  8  *  0  TO  N/2 

214  PRINT#  i?  P(B) 

215  NEXT  B 

216  ASSIGN#  1  TO  M 

217  CREATE  "GRP”,  29 

218  ASSIGN  Ml  TO  “GRP' 

219  FOR  S  »  0  TO  M 

220  PRINT#  1}  G<B) 

221  NEXT  B 

.222  ASSIGN#  1  TO  S 

225  BEEP  8  DISP  "PROGRFW  OVER” 

230  END 

231  ! 
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23S  ! 

240  !  SUBROUTINE  FKT<N,X,Y) 

250  ! 

260  !  This  rout  in*  performs  a.  FWT  of  an  input  series  in  array 
X,  The  array  Y  is  used  for  working  space. 

270  !  The  dimensions  of  X  and  Y  must  be  a  power  of  2. 

280  !  The  results  of  the  FWT  are  in  sequency  order,  positive 
290  !  phasing,  and  in  array  X.  The  subroutine  uses  a  Hadamard 
transform. 

300  ! 

310  N2  »  N/2 

320  M  ^  6  !  M  equals  the  base  2  logarithm  of  N 

330  FOR  L  *  1  TO  M 

340  N1  »  0 

350  N3  =  2^<L-D 

360  N4  »  2  X  N3 

370  N5  -  N/N4 

380  FOR  I  »  1  TO  N5 

390  N6  =  N1  ♦  i 

400  N1  »  N1  *  N3 

410  J1  »  <I-1)  X  N4 

420  J3  =  J1  +  N4  +  1 

430  FOR  J  =  N6  TO  N1 

440  J1  =  J1  +  I 

450  J2  =  J  +  N2 

460  Y<J1)  a  XCJ)  +  X<J2) 

470  J3  a  J3  -  1 

480  Y<J3)  -  X<J>-  X<J2) 

490  NEXT  J 

500  NEXT  I 

510  FOR  B  a  1  TO  N 

520  X<B)  a  Y<B) 

530  NEXT  B 
540  NEXT  L 
550  RETURN 

555  !  SUBROUTINE  PSD(N) 

560  ! 

570  ! 

580  !  This  routine  calculates  the  Walsh  Power  Spectral 
Density  Coefficients. 

590  ! 

600  P<0)  «  <F<0))^2  'first  PSD  coefficient 

610  Kal 

620  FOR  I  a  1  TO  <N/2)-l  ! coef f i ci ents  1  to  N/2-1 

630  P<I>a  <F<K))'‘2  +  <F<K+1))*2 

640  K  a  K  +  2 

650  NEXT  !I 

660  P(N/2)atF<N«l))^2  '.last  PSD  coefficient 

670  RETURN 

671  ! 

672  ! 
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680  ! 

690  !  SUBROUTINE  GROUP  SPECTRUM(M) 

700  ! 

710  !  This  subroutine  calculates  the  group  spectrum 
coeff icients- 
720  D*1 

730  G(0)  =  P<0)'^2 
740  FOR  B  =  0  TO  M-l 

750  FOR  C  »<2^B)  TO  N/2-i  STEP  <2X2''B> 

760  G<D)  *  6<0)  +  P<C)^2 

770  NEXT  !C 

780  D  =  D+1 

790  NEXT  B 

800  G<M>  =»  P<N/2)*2 

810  RETURN 

820  ! 

830  !  This  subroutine  outputs  the  Nalsh  Transform 
Coefficients. 

840  ! 

850  PRINT  -Walsh  Func t i on- ;TAe< 28) ; * Coef f i c i ent •  9  PRINT 
870  FOR  B  »  1  TO  64  ! output  the  coeffs 

880  PRINT  ;B-1  j- ,t)  "  ;TAB<20)  ;  F<B> 

900  NEXT  B 

929  RETURN 

930  ! 

940  !  Output  the  PSD  coefficients. 

958  ! 

960  PRINT 

970  PRINT  -POWER  SPECTRUM  COEFFS*  8  PRINT 

990  FOR  B  »  0  TO  N/2 

1800  PRINT  -PCs  6  •,•>=“  ;  P<8> 

1020  NEXT  B 
1040  RETURN 

1044  ! 

1045  !  Output  the  group  coefficients 

1046  ! 

1048  PRINTS  PRINT 

1049  PRINT  -GROUP  SPECTRUM  COEFFICIENTS*  8  PRINT 

1050  FOR  B  -  0  TO  M 

1060  PRINT  ■G<-jB}-)»  -iSPRINT  USING  ;GCB> 

1070  NEXT  B 
1080  RETURN 

1085  !  This  is  a  typical  DATA  line  of  numbers  representing 
TOA  TAGS 

1090  'DATA  1,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,6, 1,0, 0,0, 0,0, 0,0,1, 
0 ,0 ,0 ,0 ,0 ,0 ,0 , 1  'etc.  for  N  *  64 
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APPENDIX  B 
ORTHOGOt^LITY 

Th«  concept  of  orthogonel i ty  in  a  set  of  functions  is 
important  because  only  orthogonal  sets  of  functions  can  be 
made  to  represent  another  function  with  a  required  degree  of 
accuracy.  The  term  itself  often  brings  into  mind  the  word 
"perpendicui ar* •  and  this  thought  can  be  a  visualization  of 
the  structure  of  the  members  of  an  orthogonal  set. 

Consider  that  we  have  a  function  set*  S^^Ct),  where  n  » 

0,1, 2*3*...  The  set  is  said  to  be  orthogonal  with  weight  K 
over  the  interval  0itlT  if 

T  K  i f  n=m 

K.Sn(t)S^<  t>dt  »  (B.l) 

0  0  n^ 

with  n  and  m  being  integer  values.  If  the  constant  K  is  equal 
to  1  then  the  set  is  normalized  and  the  set  is  referred  to  as 
an  orthonormal  set. 

Mith  this  orthogonal  set,  we  may  now  represent  another 
function*  f<t>  ,  defined  over  the  interval  (0*T>*  as 

T 

f<t)  -  ZC„S„<t>  <B.2) 

8 

and  Cfi  is  a  number  that  indicates  the  value  of  the  function 
Sf,<t>.  can  be  chosen  so  as  to  mimimize  the  mean-square 

error  in  representing  f<t). 
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T 

Cj,  a  1/T  •  ’f<t)S„<t>dt  <B.3> 

e 

This  orthogonal  function  series  representation  reduces  the 
number  of  coefficients  needed  to  completely  represent  the 
signal.  C23:  pp.1-3] 


APPENDIX  C 


MODULO-2  ADDITION  AND  THE  GRAY  CODE 


Modu1o-2  addition  is  an  important  mathematical  operation 
in  Walsh  Theory.  Its  operation  is  used  in  the  definition  of 
the  Paley  ordered  Walsh  functions  with  Rademacher  functions, 
and  also  in  the  product  of  two  Walsh  functions,  namely, 

V«WL<n,t)WAL<m,t)  =WAL<n«m,t)  <C.l) 

where  n4m  indicates  modu1o-2  addition. 

Modu1o-2  addition  is  really  binary  sums  without  the  carry, 
and  obey  the  rules 

0+0=0,  0+1=1,  1+0=1,  and  1+1*0 

The  Gray  Code  is  a  binary  code  that  is  often  used  in 
communications  because  the  codes  for  successive  decimal  digits 
differs  by  only  1  bit.  It  is  not  a  weighted  code,  meaning 
that  the  decimal  value  of  a  coded  digit  cannot  be  computed  by 
a  simple  formula. 

Table  C.l.  Gray  Code  for  16  Digits. 


Decimal 

Code 

Decimal 

Code 

0 

0000 

8 

1100 

1 

0001 

9 

1101 

2 

0011 

10 

nil 

3 

00  10 

11 

1110 

4 

0  110 

12 

10  10 

5 

0111 

13 

10  1  1 

6 

0101 

14 

1001 

7 

0100 

15 

1000 
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APPB^DIX  D 


FAST  WALSH  TRANSFORM  AND  POWER  SPECTRAL  DENSITY  PLOTS 


The  -following  page-s  are  the  plots  o-f  the  FWT  and 
PSD  coe-f-f icient  o-f  the  simulated  TOA  TAG  strings  used  in  the 
thesis. 

The  plots  are  normalized  to  the  maximum  component  and  thus 
show  the  relative  value  o-f  each  component  to  this  maximum. 

The  value  o-f  the  maximum  component  is  indicated  on  the 
vertical  axis. 
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Figure  0,3.  FWT  ;  TOA  TAG  4. 
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WALSH  INDEX 
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Figure  D.7.  FHT  ; 
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Figur®  D.17.  FMT  s  TOA  TAG  22 
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WALSH  INDEX  (n) 
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Figure  D.21.  FWT  :  TOA  TAG  51. 
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WALSH  INDEX  (n) 


POWER  SPEC  DENSITY 
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Figure  D.23.  PSD  s  TOA  TAG  2. 
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SEGUENCY 


POWER  SPEC  DENSITY 
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Figure  D.25.  PSD  s  TOA  TAG  4. 
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SEQUENCY 
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PSD  :  TOA  TAG  9 
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Figure  D.43.  PSD  ;  TCA  TAG  51. 
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APPENDIX  E 


ADDITIONAL  PLOTS 

These  addi  tionjaLl  plots  are  referred  to  in  the  text 


GROUP  SPECTRUM 
TOA  TAG  5 


GROUP  NUMBER 
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Figur*  E.3.  GROUP  SPECTRUM  i  TOA  TAG  14. 
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